R Notebook
Run download_data.Rmd and
percentage_of_regional_richness.Rmd First!
library(randomForest)
library(reshape2)
library(rpart)
library(ggplot2)
library(tidyverse)
library(multcomp)
library(car)
library(VSURF)
library(boot)city_datalength(city_data$city_gdp_per_population[!is.na(city_data$city_gdp_per_population)])[1] 30length(city_data$percentage_urban_area_as_open_public_spaces[!is.na(city_data$percentage_urban_area_as_open_public_spaces)])[1] 61length(city_data$happiness_future_life[!is.na(city_data$happiness_future_life)])[1] 65length(city_data$mean_population_exposure_to_pm2_5_2019[!is.na(city_data$mean_population_exposure_to_pm2_5_2019)])[1] 131fetch_city_data_for <- function(pool_name, include_city_name = F) {
results_filename <- paste(paste('percentage_of_regional_richness__output_', pool_name, 'city', 'richness', 'intercept', sep = "_"), "csv", sep = ".")
results <- read_csv(results_filename)
joined <- left_join(city_data, results)
joined$abs_city_centre_latitude = abs(joined$city_centre_latitude)
required_columns <- c("population_growth", "rainfall_monthly_min", "rainfall_annual_average", "rainfall_monthly_max", "temperature_annual_average", "temperature_monthly_min", "temperature_monthly_max", "happiness_negative_effect", "happiness_positive_effect", "happiness_future_life", "number_of_biomes", "realm", "biome_name", "region_20km_includes_estuary", "region_50km_includes_estuary", "region_100km_includes_estuary", "city_includes_estuary", "region_20km_average_pop_density", "region_50km_average_pop_density", "region_100km_average_pop_density", "city_max_pop_density", "city_average_pop_density", "mean_population_exposure_to_pm2_5_2019", "region_20km_cultivated", "region_20km_urban", "region_50km_cultivated", "region_50km_urban", "region_100km_cultivated", "region_100km_urban", "region_20km_elevation_delta", "region_20km_mean_elevation", "region_50km_elevation_delta", "region_50km_mean_elevation", "region_100km_elevation_delta", "region_100km_mean_elevation", "city_elevation_delta", "city_mean_elevation", "urban", "shrubs", "permanent_water", "open_forest", "herbaceous_wetland", "herbaceous_vegetation", "cultivated", "closed_forest", "share_of_population_within_400m_of_open_space", "percentage_urban_area_as_streets", "percentage_urban_area_as_open_public_spaces_and_streets", "percentage_urban_area_as_open_public_spaces", "city_gdp_per_population", "city_ndvi", "city_ssm", "city_susm", "region_20km_ndvi", "region_20km_ssm", "region_20km_susm", "region_50km_ndvi", "region_50km_ssm", "region_50km_susm", "region_100km_ndvi", "region_100km_ssm", "region_100km_susm", "city_percentage_protected", "region_20km_percentage_protected", "region_50km_percentage_protected", "region_100km_percentage_protected", "city_centre_latitude", "abs_city_centre_latitude")
if (include_city_name) {
required_columns <- append(c("name"), required_columns)
}
required_columns <- append(c("response"), required_columns)
joined[,required_columns]
}| Merlin Response |
merlin_city_data <- fetch_city_data_for('merlin')
── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
cols(
name = col_character(),
response = col_double()
)
Joining, by = "name"merlin_city_datamerlin_city_data_fixed <- rfImpute(response ~ ., merlin_city_data) | Out-of-bag |
Tree | MSE %Var(y) |
300 | 21.11 117.08 |
| Out-of-bag |
Tree | MSE %Var(y) |
300 | 21.62 119.92 |
| Out-of-bag |
Tree | MSE %Var(y) |
300 | 21.99 121.97 |
| Out-of-bag |
Tree | MSE %Var(y) |
300 | 21.15 117.33 |
| Out-of-bag |
Tree | MSE %Var(y) |
300 | 21.21 117.68 |merlin_city_data_fixedggplot(merlin_city_data_fixed, aes(response)) + geom_histogram(binwidth = 2)
names(merlin_city_data_fixed) [1] "response" "population_growth"
[3] "rainfall_monthly_min" "rainfall_annual_average"
[5] "rainfall_monthly_max" "temperature_annual_average"
[7] "temperature_monthly_min" "temperature_monthly_max"
[9] "happiness_negative_effect" "happiness_positive_effect"
[11] "happiness_future_life" "number_of_biomes"
[13] "realm" "biome_name"
[15] "region_20km_includes_estuary" "region_50km_includes_estuary"
[17] "region_100km_includes_estuary" "city_includes_estuary"
[19] "region_20km_average_pop_density" "region_50km_average_pop_density"
[21] "region_100km_average_pop_density" "city_max_pop_density"
[23] "city_average_pop_density" "mean_population_exposure_to_pm2_5_2019"
[25] "region_20km_cultivated" "region_20km_urban"
[27] "region_50km_cultivated" "region_50km_urban"
[29] "region_100km_cultivated" "region_100km_urban"
[31] "region_20km_elevation_delta" "region_20km_mean_elevation"
[33] "region_50km_elevation_delta" "region_50km_mean_elevation"
[35] "region_100km_elevation_delta" "region_100km_mean_elevation"
[37] "city_elevation_delta" "city_mean_elevation"
[39] "urban" "shrubs"
[41] "permanent_water" "open_forest"
[43] "herbaceous_wetland" "herbaceous_vegetation"
[45] "cultivated" "closed_forest"
[47] "share_of_population_within_400m_of_open_space" "percentage_urban_area_as_streets"
[49] "percentage_urban_area_as_open_public_spaces_and_streets" "percentage_urban_area_as_open_public_spaces"
[51] "city_gdp_per_population" "city_ndvi"
[53] "city_ssm" "city_susm"
[55] "region_20km_ndvi" "region_20km_ssm"
[57] "region_20km_susm" "region_50km_ndvi"
[59] "region_50km_ssm" "region_50km_susm"
[61] "region_100km_ndvi" "region_100km_ssm"
[63] "region_100km_susm" "city_percentage_protected"
[65] "region_20km_percentage_protected" "region_50km_percentage_protected"
[67] "region_100km_percentage_protected" "city_centre_latitude"
[69] "abs_city_centre_latitude" merlin_response <- merlin_city_data_fixed$response
merlin_predictors <- merlin_city_data_fixed[,-1]
merlin_predictorsmerlin_interp <- VSURF(x = merlin_predictors, y = merlin_response)Thresholding step
Estimated computational time (on one core): 122.1 sec.
|
| | 0%
|
|=== | 2%
|
|====== | 4%
|
|========= | 6%
|
|=========== | 8%
|
|============== | 10%
|
|================= | 12%
|
|==================== | 14%
|
|======================= | 16%
|
|========================== | 18%
|
|============================= | 20%
|
|=============================== | 22%
|
|================================== | 24%
|
|===================================== | 26%
|
|======================================== | 28%
|
|=========================================== | 30%
|
|============================================== | 32%
|
|================================================= | 34%
|
|=================================================== | 36%
|
|====================================================== | 38%
|
|========================================================= | 40%
|
|============================================================ | 42%
|
|=============================================================== | 44%
|
|================================================================== | 46%
|
|===================================================================== | 48%
|
|======================================================================== | 50%
|
|========================================================================== | 52%
|
|============================================================================= | 54%
|
|================================================================================ | 56%
|
|=================================================================================== | 58%
|
|====================================================================================== | 60%
|
|========================================================================================= | 62%
|
|============================================================================================ | 64%
|
|============================================================================================== | 66%
|
|================================================================================================= | 68%
|
|==================================================================================================== | 70%
|
|======================================================================================================= | 72%
|
|========================================================================================================== | 74%
|
|============================================================================================================= | 76%
|
|================================================================================================================ | 78%
|
|================================================================================================================== | 80%
|
|===================================================================================================================== | 82%
|
|======================================================================================================================== | 84%
|
|=========================================================================================================================== | 86%
|
|============================================================================================================================== | 88%
|
|================================================================================================================================= | 90%
|
|==================================================================================================================================== | 92%
|
|====================================================================================================================================== | 94%
|
|========================================================================================================================================= | 96%
|
|============================================================================================================================================ | 98%
|
|===============================================================================================================================================| 100%
Interpretation step (on 46 variables)
Estimated computational time (on one core): between 58.6 sec. and 411.7 sec.
|
| | 0%
|
|=== | 2%
|
|====== | 4%
|
|========= | 7%
|
|============ | 9%
|
|================ | 11%
|
|=================== | 13%
|
|====================== | 15%
|
|========================= | 17%
|
|============================ | 20%
|
|=============================== | 22%
|
|================================== | 24%
|
|===================================== | 26%
|
|======================================== | 28%
|
|============================================ | 30%
|
|=============================================== | 33%
|
|================================================== | 35%
|
|===================================================== | 37%
|
|======================================================== | 39%
|
|=========================================================== | 41%
|
|============================================================== | 43%
|
|================================================================= | 46%
|
|==================================================================== | 48%
|
|======================================================================== | 50%
|
|=========================================================================== | 52%
|
|============================================================================== | 54%
|
|================================================================================= | 57%
|
|==================================================================================== | 59%
|
|======================================================================================= | 61%
|
|========================================================================================== | 63%
|
|============================================================================================= | 65%
|
|================================================================================================ | 67%
|
|=================================================================================================== | 70%
|
|======================================================================================================= | 72%
|
|========================================================================================================== | 74%
|
|============================================================================================================= | 76%
|
|================================================================================================================ | 78%
|
|=================================================================================================================== | 80%
|
|====================================================================================================================== | 83%
|
|========================================================================================================================= | 85%
|
|============================================================================================================================ | 87%
|
|=============================================================================================================================== | 89%
|
|=================================================================================================================================== | 91%
|
|====================================================================================================================================== | 93%
|
|========================================================================================================================================= | 96%
|
|============================================================================================================================================ | 98%
|
|===============================================================================================================================================| 100%
Prediction step (on 3 variables)
Maximum estimated computational time (on one core): 3.6 sec.
|
| | 0%
|
|================================================ | 33%
|
|=============================================================================================== | 67%
|
|===============================================================================================================================================| 100%names(merlin_predictors[,merlin_interp$varselect.interp])[1] "abs_city_centre_latitude" "region_50km_elevation_delta" "biome_name" | Birdlife Response |
birdlife_city_data <- fetch_city_data_for('birdlife')
── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
cols(
name = col_character(),
response = col_double()
)
Joining, by = "name"birdlife_city_dataggplot(birdlife_city_data, aes(response)) + geom_histogram(binwidth = 1)
birdlife_city_data_fixed <- rfImpute(response ~ ., birdlife_city_data) | Out-of-bag |
Tree | MSE %Var(y) |
300 | 5.621 88.99 |
| Out-of-bag |
Tree | MSE %Var(y) |
300 | 5.76 91.18 |
| Out-of-bag |
Tree | MSE %Var(y) |
300 | 5.793 91.71 |
| Out-of-bag |
Tree | MSE %Var(y) |
300 | 5.846 92.54 |
| Out-of-bag |
Tree | MSE %Var(y) |
300 | 5.701 90.24 |birdlife_city_data_fixednames(birdlife_city_data_fixed) [1] "response" "population_growth"
[3] "rainfall_monthly_min" "rainfall_annual_average"
[5] "rainfall_monthly_max" "temperature_annual_average"
[7] "temperature_monthly_min" "temperature_monthly_max"
[9] "happiness_negative_effect" "happiness_positive_effect"
[11] "happiness_future_life" "number_of_biomes"
[13] "realm" "biome_name"
[15] "region_20km_includes_estuary" "region_50km_includes_estuary"
[17] "region_100km_includes_estuary" "city_includes_estuary"
[19] "region_20km_average_pop_density" "region_50km_average_pop_density"
[21] "region_100km_average_pop_density" "city_max_pop_density"
[23] "city_average_pop_density" "mean_population_exposure_to_pm2_5_2019"
[25] "region_20km_cultivated" "region_20km_urban"
[27] "region_50km_cultivated" "region_50km_urban"
[29] "region_100km_cultivated" "region_100km_urban"
[31] "region_20km_elevation_delta" "region_20km_mean_elevation"
[33] "region_50km_elevation_delta" "region_50km_mean_elevation"
[35] "region_100km_elevation_delta" "region_100km_mean_elevation"
[37] "city_elevation_delta" "city_mean_elevation"
[39] "urban" "shrubs"
[41] "permanent_water" "open_forest"
[43] "herbaceous_wetland" "herbaceous_vegetation"
[45] "cultivated" "closed_forest"
[47] "share_of_population_within_400m_of_open_space" "percentage_urban_area_as_streets"
[49] "percentage_urban_area_as_open_public_spaces_and_streets" "percentage_urban_area_as_open_public_spaces"
[51] "city_gdp_per_population" "city_ndvi"
[53] "city_ssm" "city_susm"
[55] "region_20km_ndvi" "region_20km_ssm"
[57] "region_20km_susm" "region_50km_ndvi"
[59] "region_50km_ssm" "region_50km_susm"
[61] "region_100km_ndvi" "region_100km_ssm"
[63] "region_100km_susm" "city_percentage_protected"
[65] "region_20km_percentage_protected" "region_50km_percentage_protected"
[67] "region_100km_percentage_protected" "city_centre_latitude"
[69] "abs_city_centre_latitude" birdlife_response <- birdlife_city_data_fixed$response
birdlife_predictors <- birdlife_city_data_fixed[,-1]
birdlife_predictorsbirdlife_interp <- VSURF(x = birdlife_predictors, y = birdlife_response)Thresholding step
Estimated computational time (on one core): 106.5 sec.
|
| | 0%
|
|=== | 2%
|
|====== | 4%
|
|========= | 6%
|
|=========== | 8%
|
|============== | 10%
|
|================= | 12%
|
|==================== | 14%
|
|======================= | 16%
|
|========================== | 18%
|
|============================= | 20%
|
|=============================== | 22%
|
|================================== | 24%
|
|===================================== | 26%
|
|======================================== | 28%
|
|=========================================== | 30%
|
|============================================== | 32%
|
|================================================= | 34%
|
|=================================================== | 36%
|
|====================================================== | 38%
|
|========================================================= | 40%
|
|============================================================ | 42%
|
|=============================================================== | 44%
|
|================================================================== | 46%
|
|===================================================================== | 48%
|
|======================================================================== | 50%
|
|========================================================================== | 52%
|
|============================================================================= | 54%
|
|================================================================================ | 56%
|
|=================================================================================== | 58%
|
|====================================================================================== | 60%
|
|========================================================================================= | 62%
|
|============================================================================================ | 64%
|
|============================================================================================== | 66%
|
|================================================================================================= | 68%
|
|==================================================================================================== | 70%
|
|======================================================================================================= | 72%
|
|========================================================================================================== | 74%
|
|============================================================================================================= | 76%
|
|================================================================================================================ | 78%
|
|================================================================================================================== | 80%
|
|===================================================================================================================== | 82%
|
|======================================================================================================================== | 84%
|
|=========================================================================================================================== | 86%
|
|============================================================================================================================== | 88%
|
|================================================================================================================================= | 90%
|
|==================================================================================================================================== | 92%
|
|====================================================================================================================================== | 94%
|
|========================================================================================================================================= | 96%
|
|============================================================================================================================================ | 98%
|
|===============================================================================================================================================| 100%
Interpretation step (on 55 variables)
Estimated computational time (on one core): between 72.9 sec. and 534.9 sec.
|
| | 0%
|
|=== | 2%
|
|===== | 4%
|
|======== | 5%
|
|========== | 7%
|
|============= | 9%
|
|================ | 11%
|
|================== | 13%
|
|===================== | 15%
|
|======================= | 16%
|
|========================== | 18%
|
|============================= | 20%
|
|=============================== | 22%
|
|================================== | 24%
|
|==================================== | 25%
|
|======================================= | 27%
|
|========================================== | 29%
|
|============================================ | 31%
|
|=============================================== | 33%
|
|================================================= | 35%
|
|==================================================== | 36%
|
|======================================================= | 38%
|
|========================================================= | 40%
|
|============================================================ | 42%
|
|============================================================== | 44%
|
|================================================================= | 45%
|
|==================================================================== | 47%
|
|====================================================================== | 49%
|
|========================================================================= | 51%
|
|=========================================================================== | 53%
|
|============================================================================== | 55%
|
|================================================================================= | 56%
|
|=================================================================================== | 58%
|
|====================================================================================== | 60%
|
|======================================================================================== | 62%
|
|=========================================================================================== | 64%
|
|============================================================================================== | 65%
|
|================================================================================================ | 67%
|
|=================================================================================================== | 69%
|
|===================================================================================================== | 71%
|
|======================================================================================================== | 73%
|
|=========================================================================================================== | 75%
|
|============================================================================================================= | 76%
|
|================================================================================================================ | 78%
|
|================================================================================================================== | 80%
|
|===================================================================================================================== | 82%
|
|======================================================================================================================== | 84%
|
|========================================================================================================================== | 85%
|
|============================================================================================================================= | 87%
|
|=============================================================================================================================== | 89%
|
|================================================================================================================================== | 91%
|
|===================================================================================================================================== | 93%
|
|======================================================================================================================================= | 95%
|
|========================================================================================================================================== | 96%
|
|============================================================================================================================================ | 98%
|
|===============================================================================================================================================| 100%
Prediction step (on 2 variables)
Maximum estimated computational time (on one core): 2.6 sec.
|
| | 0%
|
|======================================================================== | 50%
|
|===============================================================================================================================================| 100%names(birdlife_predictors[,birdlife_interp$varselect.interp])[1] "population_growth" "region_50km_ssm" | So…. |
|---|
| Merlin: “abs_city_centre_latitude” “region_50km_elevation_delta” “biome_name” “region_50km_ssm” [5] “region_100km_ssm” “region_20km_elevation_delta” “region_20km_urban” “region_100km_elevation_delta” [9] “temperature_annual_average” “city_ndvi” “city_gdp_per_population” “shrubs” [13] “permanent_water” “city_centre_latitude” “region_20km_cultivated” Birdlife: “population_growth” “region_50km_ssm” |
Try Modelling
merlin_city_data_named <- fetch_city_data_for('merlin', T)
── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
cols(
name = col_character(),
response = col_double()
)
Joining, by = "name"birdlife_city_data_named <- fetch_city_data_for('birdlife', T)
── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
cols(
name = col_character(),
response = col_double()
)
Joining, by = "name"| Use cross validation and dropping terms to find best model |
full model: response ~ region_100km_ssm + region_20km_elevation_delta + region_20km_urban + region_100km_elevation_delta + temperature_annual_average + city_ndvi + city_gdp_per_population + shrubs + permanent_water + city_centre_latitude + region_20km_cultivated
Merlin data set
merlin_city_data_fixed_no_boreal <- merlin_city_data_fixed[merlin_city_data_fixed$biome_name != 'Boreal Forests/Taiga',]cv.glm(merlin_city_data_fixed_no_boreal, glm(data = merlin_city_data_fixed_no_boreal, formula = response ~ region_100km_ssm + region_20km_elevation_delta + region_20km_urban + region_100km_elevation_delta + temperature_annual_average + city_ndvi + city_gdp_per_population + shrubs + permanent_water + city_centre_latitude + region_20km_cultivated))$delta[1][1] 19.88982– CVE 19.72841 – Can we drop one?
cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_20km_elevation_delta + region_20km_urban + region_100km_elevation_delta + temperature_annual_average + city_ndvi + city_gdp_per_population + shrubs + permanent_water + city_centre_latitude + region_20km_cultivated, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 19.97547cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_urban + region_100km_elevation_delta + temperature_annual_average + city_ndvi + city_gdp_per_population + shrubs + permanent_water + city_centre_latitude + region_20km_cultivated, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 19.54162cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_elevation_delta + region_100km_elevation_delta + temperature_annual_average + city_ndvi + city_gdp_per_population + shrubs + permanent_water + city_centre_latitude + region_20km_cultivated, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 19.77657cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_elevation_delta + region_20km_urban + temperature_annual_average + city_ndvi + city_gdp_per_population + shrubs + permanent_water + city_centre_latitude + region_20km_cultivated, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 19.83879cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_elevation_delta + region_20km_urban + region_100km_elevation_delta + city_ndvi + city_gdp_per_population + shrubs + permanent_water + city_centre_latitude + region_20km_cultivated, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 19.63049cv.glm(merlin_city_data_fixed_no_boreal, glm(data = merlin_city_data_fixed_no_boreal, formula = response ~ region_100km_ssm + region_20km_elevation_delta + region_20km_urban + region_100km_elevation_delta + temperature_annual_average + city_gdp_per_population + shrubs + permanent_water + city_centre_latitude + region_20km_cultivated))$delta[1][1] 19.44958cv.glm(merlin_city_data_fixed_no_boreal, glm(data = merlin_city_data_fixed_no_boreal, formula = response ~ region_100km_ssm + region_20km_elevation_delta + region_20km_urban + region_100km_elevation_delta + temperature_annual_average + city_ndvi + shrubs + permanent_water + city_centre_latitude + region_20km_cultivated))$delta[1][1] 19.74441cv.glm(merlin_city_data_fixed_no_boreal, glm(data = merlin_city_data_fixed_no_boreal, formula = response ~ region_100km_ssm + region_20km_elevation_delta + region_20km_urban + region_100km_elevation_delta + temperature_annual_average + city_ndvi + city_gdp_per_population + permanent_water + city_centre_latitude + region_20km_cultivated))$delta[1][1] 19.75178cv.glm(merlin_city_data_fixed_no_boreal, glm(data = merlin_city_data_fixed_no_boreal, formula = response ~ region_100km_ssm + region_20km_elevation_delta + region_20km_urban + region_100km_elevation_delta + temperature_annual_average + city_ndvi + city_gdp_per_population + shrubs + city_centre_latitude + region_20km_cultivated))$delta[1][1] 19.70767cv.glm(merlin_city_data_fixed_no_boreal, glm(data = merlin_city_data_fixed_no_boreal, formula = response ~ region_100km_ssm + region_20km_elevation_delta + region_20km_urban + region_100km_elevation_delta + temperature_annual_average + city_ndvi + city_gdp_per_population + shrubs + permanent_water + region_20km_cultivated))$delta[1][1] 19.58658cv.glm(merlin_city_data_fixed_no_boreal, glm(data = merlin_city_data_fixed_no_boreal, formula = response ~ region_100km_ssm + region_20km_elevation_delta + region_20km_urban + region_100km_elevation_delta + temperature_annual_average + city_ndvi + city_gdp_per_population + shrubs + permanent_water + city_centre_latitude))$delta[1][1] 19.62306– drop city_ndvi to give smaller CVE of 19.35 – can we drop another?
cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_20km_elevation_delta + region_20km_urban + region_100km_elevation_delta + temperature_annual_average + city_gdp_per_population + shrubs + permanent_water + city_centre_latitude + region_20km_cultivated, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 19.58667cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_urban + region_100km_elevation_delta + temperature_annual_average + city_gdp_per_population + shrubs + permanent_water + city_centre_latitude + region_20km_cultivated, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 19.12746cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_elevation_delta + region_100km_elevation_delta + temperature_annual_average + city_gdp_per_population + shrubs + permanent_water + city_centre_latitude + region_20km_cultivated, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 19.35618cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_elevation_delta + region_20km_urban + region_100km_elevation_delta + city_gdp_per_population + shrubs + permanent_water + city_centre_latitude + region_20km_cultivated, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 19.20135cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_elevation_delta + region_20km_urban + region_100km_elevation_delta + temperature_annual_average + shrubs + permanent_water + city_centre_latitude + region_20km_cultivated, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 19.31519cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_elevation_delta + region_20km_urban + region_100km_elevation_delta + temperature_annual_average + city_gdp_per_population + permanent_water + city_centre_latitude + region_20km_cultivated, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 19.32433cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_elevation_delta + region_20km_urban + region_100km_elevation_delta + temperature_annual_average + city_gdp_per_population + shrubs + city_centre_latitude + region_20km_cultivated, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 19.28881cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_elevation_delta + region_20km_urban + region_100km_elevation_delta + temperature_annual_average + city_gdp_per_population + shrubs + permanent_water + region_20km_cultivated, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 19.1589cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_elevation_delta + region_20km_urban + region_100km_elevation_delta + temperature_annual_average + city_gdp_per_population + shrubs + permanent_water + city_centre_latitude, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 19.18252– drop city_centre_latitude to give smaller CVE of 19.06 – can we drop another?
cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_20km_elevation_delta + region_20km_urban + region_100km_elevation_delta + temperature_annual_average + city_gdp_per_population + shrubs + permanent_water + region_20km_cultivated, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 19.31336cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_urban + region_100km_elevation_delta + temperature_annual_average + city_gdp_per_population + shrubs + permanent_water + region_20km_cultivated, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.85153cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_elevation_delta + region_100km_elevation_delta + temperature_annual_average + city_gdp_per_population + shrubs + permanent_water + region_20km_cultivated, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 19.07111cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_elevation_delta + region_20km_urban + temperature_annual_average + city_gdp_per_population + shrubs + permanent_water + region_20km_cultivated, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 19.12285cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_elevation_delta + region_20km_urban + region_100km_elevation_delta + city_gdp_per_population + shrubs + permanent_water + region_20km_cultivated, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.93829cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_elevation_delta + region_20km_urban + region_100km_elevation_delta + temperature_annual_average + shrubs + permanent_water + region_20km_cultivated, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 19.02621cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_elevation_delta + region_20km_urban + region_100km_elevation_delta + temperature_annual_average + city_gdp_per_population + permanent_water + region_20km_cultivated, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 19.04182cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_elevation_delta + region_20km_urban + region_100km_elevation_delta + temperature_annual_average + city_gdp_per_population + shrubs + region_20km_cultivated, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 19.00628cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_elevation_delta + region_20km_urban + region_100km_elevation_delta + temperature_annual_average + city_gdp_per_population + shrubs + permanent_water, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.91404– drop region_20km_elevation_delta to give smaller CVE of 18.76 – can we drop another?
cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_20km_urban + region_100km_elevation_delta + temperature_annual_average + city_gdp_per_population + shrubs + permanent_water + region_20km_cultivated, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.97297cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_100km_elevation_delta + temperature_annual_average + city_gdp_per_population + shrubs + permanent_water + region_20km_cultivated, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.74933cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_urban + temperature_annual_average + city_gdp_per_population + shrubs + permanent_water + region_20km_cultivated, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.89437cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_urban + region_100km_elevation_delta + city_gdp_per_population + shrubs + permanent_water + region_20km_cultivated, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.63393cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_urban + region_100km_elevation_delta + temperature_annual_average + shrubs + permanent_water + region_20km_cultivated, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.72148cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_urban + region_100km_elevation_delta + temperature_annual_average + city_gdp_per_population + permanent_water + region_20km_cultivated, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.73854cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_urban + region_100km_elevation_delta + temperature_annual_average + city_gdp_per_population + shrubs + region_20km_cultivated, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.69534cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_urban + region_100km_elevation_delta + temperature_annual_average + city_gdp_per_population + shrubs + permanent_water, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.61894– drop region_20km_cultivated to give smaller CVE of 18.54 – can we drop another?
cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_20km_urban + region_100km_elevation_delta + temperature_annual_average + city_gdp_per_population + shrubs + permanent_water, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.79968cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_100km_elevation_delta + temperature_annual_average + city_gdp_per_population + shrubs + permanent_water, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.50997cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_urban + temperature_annual_average + city_gdp_per_population + shrubs + permanent_water, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.68066cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_urban + region_100km_elevation_delta + city_gdp_per_population + shrubs + permanent_water, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.40791cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_urban + region_100km_elevation_delta + temperature_annual_average + shrubs + permanent_water, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.48533cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_urban + region_100km_elevation_delta + temperature_annual_average + city_gdp_per_population + permanent_water, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.53006cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_urban + region_100km_elevation_delta + temperature_annual_average + city_gdp_per_population + shrubs, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.46614– drop permanent_water to give smaller CVE of 18.34 – can we drop another?
cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_20km_urban + region_100km_elevation_delta + temperature_annual_average + city_gdp_per_population + shrubs, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.57444cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_100km_elevation_delta + temperature_annual_average + city_gdp_per_population + shrubs, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.37669cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_urban + temperature_annual_average + city_gdp_per_population + shrubs, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.58396cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_urban + region_100km_elevation_delta + city_gdp_per_population + shrubs, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.26614cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_urban + region_100km_elevation_delta + temperature_annual_average + shrubs, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.33385cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_urban + region_100km_elevation_delta + temperature_annual_average + city_gdp_per_population, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.39355– drop temperature_annual_average to give smaller CVE of 18.14 – can we drop another?
cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_20km_urban + region_100km_elevation_delta + city_gdp_per_population + shrubs, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.37121cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_100km_elevation_delta + city_gdp_per_population + shrubs, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.13177cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_urban + city_gdp_per_population + shrubs, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.36207cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_urban + region_100km_elevation_delta + shrubs, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.14642cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_20km_urban + region_100km_elevation_delta + city_gdp_per_population, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.1946– drop region_20km_urban to give smaller CVE of 18.03 – can we drop another?
cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_elevation_delta + city_gdp_per_population + shrubs, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.33203cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + city_gdp_per_population + shrubs, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.22241cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_100km_elevation_delta + shrubs, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.01086cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_100km_elevation_delta + city_gdp_per_population, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.05821– drop shrubs to give smaller CVE of 17.95 – can we drop another?
cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_elevation_delta + city_gdp_per_population, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.22109cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + city_gdp_per_population, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.25476cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_100km_elevation_delta, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 17.94321– drop city_gdp_per_population to give smaller CVE of 17.94 – can we drop another?
cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_elevation_delta, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.06258cv.glm(merlin_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm, data = merlin_city_data_fixed_no_boreal))$delta[1][1] 18.19515– No
| – best model with region_100km_ssm + region_100km_elevation_delta (CV error 17.91216) |
summary(glm(data = merlin_city_data_fixed, formula = response ~ region_100km_ssm + region_100km_elevation_delta))
Call:
glm(formula = response ~ region_100km_ssm + region_100km_elevation_delta,
data = merlin_city_data_fixed)
Deviance Residuals:
Min 1Q Median 3Q Max
-7.6625 -3.0246 -0.4496 1.9868 16.9640
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 2.6877517 1.1338124 2.371 0.0192 *
region_100km_ssm -0.1327133 0.0698210 -1.901 0.0595 .
region_100km_elevation_delta -0.0005261 0.0002944 -1.787 0.0762 .
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 17.46142)
Null deviance: 2469.6 on 136 degrees of freedom
Residual deviance: 2339.8 on 134 degrees of freedom
AIC: 785.57
Number of Fisher Scoring iterations: 2reg_merlin = glm(data = merlin_city_data_fixed, formula = response ~ region_100km_ssm + region_100km_elevation_delta)
with(summary(reg_merlin), 1 - deviance/null.deviance)[1] 0.05255Birdlife data set
birdlife_city_data_fixed_no_boreal <- birdlife_city_data_fixed[birdlife_city_data_fixed$biome_name != 'Boreal Forests/Taiga',]cv.glm(birdlife_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_50km_ssm + region_50km_elevation_delta + biome_name + population_growth, data = birdlife_city_data_fixed_no_boreal))$delta[1][1] 6.899862– can we drop a variable?
cv.glm(birdlife_city_data_fixed_no_boreal, glm(formula = response ~ region_50km_ssm + region_50km_elevation_delta + biome_name + population_growth, data = birdlife_city_data_fixed_no_boreal))$delta[1][1] 6.768164cv.glm(birdlife_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_50km_elevation_delta + biome_name + population_growth, data = birdlife_city_data_fixed_no_boreal))$delta[1][1] 6.752211cv.glm(birdlife_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_50km_ssm + biome_name + population_growth, data = birdlife_city_data_fixed_no_boreal))$delta[1][1] 6.989636cv.glm(birdlife_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_50km_ssm + region_50km_elevation_delta + population_growth, data = birdlife_city_data_fixed_no_boreal))$delta[1][1] 6.503421cv.glm(birdlife_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_50km_ssm + region_50km_elevation_delta + biome_name, data = birdlife_city_data_fixed_no_boreal))$delta[1][1] 6.780578– drop biome_name to give CVE of 6.503421 – can we drop another?
cv.glm(birdlife_city_data_fixed_no_boreal, glm(formula = response ~ region_50km_ssm + region_50km_elevation_delta + population_growth, data = birdlife_city_data_fixed_no_boreal))$delta[1][1] 6.417311cv.glm(birdlife_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_50km_elevation_delta + population_growth, data = birdlife_city_data_fixed_no_boreal))$delta[1][1] 6.426562cv.glm(birdlife_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_50km_ssm + population_growth, data = birdlife_city_data_fixed_no_boreal))$delta[1][1] 6.430742cv.glm(birdlife_city_data_fixed_no_boreal, glm(formula = response ~ region_100km_ssm + region_50km_ssm + region_50km_elevation_delta, data = birdlife_city_data_fixed_no_boreal))$delta[1][1] 6.439714– drop region_100km_ssm to give CVE of 6.417311 – can we drop another?
cv.glm(birdlife_city_data_fixed_no_boreal, glm(formula = response ~ region_50km_elevation_delta + population_growth, data = birdlife_city_data_fixed_no_boreal))$delta[1][1] 6.535285cv.glm(birdlife_city_data_fixed_no_boreal, glm(formula = response ~ region_50km_ssm + population_growth, data = birdlife_city_data_fixed_no_boreal))$delta[1][1] 6.342025cv.glm(birdlife_city_data_fixed_no_boreal, glm(formula = response ~ region_50km_ssm + region_50km_elevation_delta, data = birdlife_city_data_fixed_no_boreal))$delta[1][1] 6.352664– drop region_50km_elevation_delta to give CVE of 6.342025 – can we drop another?
cv.glm(birdlife_city_data_fixed_no_boreal, glm(formula = response ~ population_growth, data = birdlife_city_data_fixed_no_boreal))$delta[1][1] 6.464699cv.glm(birdlife_city_data_fixed_no_boreal, glm(formula = response ~ region_50km_ssm, data = birdlife_city_data_fixed_no_boreal))$delta[1][1] 6.291299– drop city_gdp_per_population to give CVE of 6.291299 – is this better than no variable?
cv.glm(birdlife_city_data_fixed_no_boreal, glm(formula = response ~ 1, data = birdlife_city_data_fixed_no_boreal))$delta[1][1] 6.395701– yes, just!
| – so best model with birdlife is region_50km_ssm |
summary(glm(data = birdlife_city_data_fixed, formula = response ~ region_50km_ssm))
Call:
glm(formula = response ~ region_50km_ssm, data = birdlife_city_data_fixed)
Deviance Residuals:
Min 1Q Median 3Q Max
-4.5353 -1.5461 -0.4124 1.3071 10.7572
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.26916 0.65041 1.951 0.0531 .
region_50km_ssm -0.08499 0.04115 -2.065 0.0408 *
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 6.214378)
Null deviance: 865.45 on 136 degrees of freedom
Residual deviance: 838.94 on 135 degrees of freedom
AIC: 643.06
Number of Fisher Scoring iterations: 2reg_birdlife = glm(data = birdlife_city_data_fixed, formula = response ~ region_50km_ssm)
with(summary(reg_birdlife), 1 - deviance/null.deviance)[1] 0.03062471ggplot(birdlife_city_data_named, aes(x = region_50km_ssm, y = response)) + geom_point() + geom_smooth(method = "glm")`geom_smooth()` using formula 'y ~ x'
| Check birdlife model fit |
birdlife.fit <- glm(data = birdlife_city_data_fixed, formula = response ~ region_50km_ssm)
summary(birdlife.fit)
Call:
glm(formula = response ~ region_50km_ssm, data = birdlife_city_data_fixed)
Deviance Residuals:
Min 1Q Median 3Q Max
-4.5353 -1.5461 -0.4124 1.3071 10.7572
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.26916 0.65041 1.951 0.0531 .
region_50km_ssm -0.08499 0.04115 -2.065 0.0408 *
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 6.214378)
Null deviance: 865.45 on 136 degrees of freedom
Residual deviance: 838.94 on 135 degrees of freedom
AIC: 643.06
Number of Fisher Scoring iterations: 2with(summary(birdlife.fit), 1 - deviance/null.deviance)[1] 0.03062471plot(birdlife.fit)
birdlife_city_data_fixed_no_boreal[c(16, 53, 72), c("region_50km_ssm")][1] 18.451180 9.961682 11.644862city_data[c(16, 53, 72), c("name", "region_50km_ssm")]dat <- predict(glm(formula = response ~ region_50km_ssm, data = birdlife_city_data_named), se.fit=T)
outside_se <- birdlife_city_data_named[birdlife_city_data_named$response < dat$fit - 15* dat$se.fit | birdlife_city_data_named$response > dat$fit + 15 * dat$se.fit,]
ggplot(birdlife_city_data_named, aes(x = region_50km_ssm, y = response)) +
geom_point(size=1) +
geom_smooth(method = "glm") +
geom_text(aes(label = name), data = outside_se, size = 3, position = "dodge", vjust = "inward", hjust = "inward", color = "red", angle=-15) +
geom_point(data = outside_se, color = "red") +
theme_bw() +
ylab("City Random Effect Intercept") + xlab("Regional (50km) SSM") + labs(title = "Birdlife")`geom_smooth()` using formula 'y ~ x'
Warning: Width not defined. Set with `position_dodge(width = ?)`ggsave("city_effect_richness__output__birdlife.jpg")Saving 7.29 x 4.51 in image
`geom_smooth()` using formula 'y ~ x'
Warning: Width not defined. Set with `position_dodge(width = ?)`
How much variation have we explained?
birdlife_city_data_named$residuals <- resid(birdlife.fit)
ggplot(birdlife_city_data_named, aes(y = response, x = residuals)) +
geom_point() +
geom_smooth(method = "lm", se = F) +
labs(title = "Birdlife", subtitle = paste("Correlation", cor(birdlife_city_data_named$residuals, birdlife_city_data_named$response))) +
theme_bw()`geom_smooth()` using formula 'y ~ x'
| Check Merlin model fit |
merlin.fit <- glm(data = merlin_city_data_named, formula = response ~ region_100km_ssm + region_50km_elevation_delta)
summary(merlin.fit)
Call:
glm(formula = response ~ region_100km_ssm + region_50km_elevation_delta,
data = merlin_city_data_named)
Deviance Residuals:
Min 1Q Median 3Q Max
-7.6599 -2.9987 -0.5524 1.7449 16.9143
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 2.6809304 1.1210300 2.391 0.0182 *
region_100km_ssm -0.1331207 0.0695604 -1.914 0.0578 .
region_50km_elevation_delta -0.0006899 0.0003461 -1.994 0.0482 *
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 17.36262)
Null deviance: 2469.6 on 136 degrees of freedom
Residual deviance: 2326.6 on 134 degrees of freedom
AIC: 784.8
Number of Fisher Scoring iterations: 2with(summary(merlin.fit), 1 - deviance/null.deviance)[1] 0.05791102plot(merlin.fit)
merlin_city_data_fixed_no_boreal[c(24, 42, 72), c("region_100km_ssm", "region_50km_elevation_delta")]city_data[c(24, 42, 72), c("name", "region_100km_ssm", "region_50km_elevation_delta")]ggplot(merlin_city_data_named, aes(x = region_100km_ssm, y = response)) +
geom_point(aes(size = region_50km_elevation_delta)) +
geom_smooth(method = "glm") +
theme_bw() +
theme(legend.position="bottom") +
ylab("City Random Effect Intercept") + xlab("Regional (100km) SSM") + labs(title = "eBird") + guides(size=guide_legend(title="Regional (50km) Elevation Delta"))`geom_smooth()` using formula 'y ~ x'ggsave("city_effect_richness__output__merlin.jpg")Saving 7.29 x 4.51 in image
`geom_smooth()` using formula 'y ~ x'
How much variation have we explained?
merlin_city_data_named$residuals <- resid(merlin.fit)
ggplot(merlin_city_data_named, aes(y = response, x = residuals)) +
geom_point() +
geom_smooth(method = "lm", se = F) +
labs(title = "Merlin", subtitle = paste("Correlation", cor(merlin_city_data_named$residuals, merlin_city_data_named$response))) +
theme_bw()`geom_smooth()` using formula 'y ~ x'
| Check AIC |
AIC(
glm(data = merlin_city_data_fixed, formula = response ~ region_100km_ssm + region_50km_ssm + region_50km_elevation_delta + biome_name + population_growth),
glm(data = merlin_city_data_fixed, formula = response ~ region_100km_ssm + region_50km_elevation_delta)
)AIC(
glm(data = birdlife_city_data_fixed, formula = response ~ region_100km_ssm + region_50km_ssm + region_50km_elevation_delta + biome_name + population_growth),
glm(data = birdlife_city_data_fixed, formula = response ~ region_50km_ssm)
)| LDG |
birdlife_data_with_lat = fetch_city_data_for('birdlife', include_city_name = T)
── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
cols(
name = col_character(),
response = col_double()
)
Joining, by = "name"model2 <- glm(formula = response ~ I(city_centre_latitude^2), data = birdlife_data_with_lat)
dat2 <- predict(model2, se.fit=T)
summary(model2)
Call:
glm(formula = response ~ I(city_centre_latitude^2), data = birdlife_data_with_lat)
Deviance Residuals:
Min 1Q Median 3Q Max
-4.7461 -1.4878 -0.4638 1.2399 9.5839
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.4426687 0.3308024 -1.338 0.1831
I(city_centre_latitude^2) 0.0004780 0.0002725 1.754 0.0817 .
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 6.267835)
Null deviance: 865.45 on 136 degrees of freedom
Residual deviance: 846.16 on 135 degrees of freedom
AIC: 644.23
Number of Fisher Scoring iterations: 2with(summary(model2), 1 - deviance/null.deviance)[1] 0.02228615outside_se <- birdlife_data_with_lat[birdlife_data_with_lat$response < dat2$fit - 15* dat2$se.fit | birdlife_data_with_lat$response > dat2$fit + 15 * dat2$se.fit,]
ggplot(birdlife_data_with_lat, aes(x = city_centre_latitude, y = response)) +
geom_point(size=1) +
geom_smooth(method = "glm", formula = y ~ I(x^2)) +
geom_text(aes(label = name), data = outside_se, size = 3, position = "dodge", vjust = "inward", hjust = "inward", color = "red", angle=-15) +
geom_point(data = outside_se, color = "red") +
theme_bw() +
ylab("City Random Effect Intercept") + xlab("Latitude (of city centre)") + labs(title = "Birdlife")Warning: Width not defined. Set with `position_dodge(width = ?)`ggsave('city_effect_richness__output__ldg_birdlife.jpg')Saving 7.29 x 4.51 in image
Warning: Width not defined. Set with `position_dodge(width = ?)`
merlin_data_with_lat = fetch_city_data_for('merlin', include_city_name = T)
── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
cols(
name = col_character(),
response = col_double()
)
Joining, by = "name"model2 <- glm(formula = response ~ I(city_centre_latitude^2), data = merlin_data_with_lat)
dat2 <- predict(model2, se.fit=T)
summary(model2)
Call:
glm(formula = response ~ I(city_centre_latitude^2), data = merlin_data_with_lat)
Deviance Residuals:
Min 1Q Median 3Q Max
-9.209 -2.825 -0.388 1.463 18.438
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.201e-02 5.651e-01 0.021 0.983
I(city_centre_latitude^2) -1.296e-05 4.655e-04 -0.028 0.978
(Dispersion parameter for gaussian family taken to be 18.29329)
Null deviance: 2469.6 on 136 degrees of freedom
Residual deviance: 2469.6 on 135 degrees of freedom
AIC: 790.97
Number of Fisher Scoring iterations: 2with(summary(model2), 1 - deviance/null.deviance)[1] 5.74457e-06outside_se <- merlin_data_with_lat[merlin_data_with_lat$response < dat2$fit - 15* dat2$se.fit | merlin_data_with_lat$response > dat2$fit + 15 * dat2$se.fit,]
ggplot(merlin_data_with_lat, aes(x = city_centre_latitude, y = response)) +
geom_point(size=1) +
geom_smooth(method = "glm", formula = y ~ I(x^2)) +
geom_text(aes(label = name), data = outside_se, size = 3, position = "dodge", vjust = "inward", hjust = "inward", color = "red", angle=-15) +
geom_point(data = outside_se, color = "red") +
theme_bw() +
ylab("City Random Effect Intercept") + xlab("Latitude (of city centre)") + labs(title = "eBird")Warning: Width not defined. Set with `position_dodge(width = ?)`ggsave('city_effect_richness__output__ldg_merlin.jpg')Saving 7.29 x 4.51 in image
Warning: Width not defined. Set with `position_dodge(width = ?)`
LS0tCnRpdGxlOiAiUiBOb3RlYm9vayIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQpSdW4gYGRvd25sb2FkX2RhdGEuUm1kYCBhbmQgYHBlcmNlbnRhZ2Vfb2ZfcmVnaW9uYWxfcmljaG5lc3MuUm1kYCBGaXJzdCEKCmBgYHtyIHNldHVwfQpsaWJyYXJ5KHJhbmRvbUZvcmVzdCkKbGlicmFyeShyZXNoYXBlMikKbGlicmFyeShycGFydCkKbGlicmFyeShnZ3Bsb3QyKQpsaWJyYXJ5KHRpZHl2ZXJzZSkKCmxpYnJhcnkobXVsdGNvbXApCmxpYnJhcnkoY2FyKQpsaWJyYXJ5KFZTVVJGKQoKbGlicmFyeShib290KQpgYGAKCmBgYHtyfQpjaXR5X2RhdGEKYGBgCgpgYGB7cn0KbGVuZ3RoKGNpdHlfZGF0YSRjaXR5X2dkcF9wZXJfcG9wdWxhdGlvblshaXMubmEoY2l0eV9kYXRhJGNpdHlfZ2RwX3Blcl9wb3B1bGF0aW9uKV0pCmxlbmd0aChjaXR5X2RhdGEkcGVyY2VudGFnZV91cmJhbl9hcmVhX2FzX29wZW5fcHVibGljX3NwYWNlc1shaXMubmEoY2l0eV9kYXRhJHBlcmNlbnRhZ2VfdXJiYW5fYXJlYV9hc19vcGVuX3B1YmxpY19zcGFjZXMpXSkKbGVuZ3RoKGNpdHlfZGF0YSRoYXBwaW5lc3NfZnV0dXJlX2xpZmVbIWlzLm5hKGNpdHlfZGF0YSRoYXBwaW5lc3NfZnV0dXJlX2xpZmUpXSkKbGVuZ3RoKGNpdHlfZGF0YSRtZWFuX3BvcHVsYXRpb25fZXhwb3N1cmVfdG9fcG0yXzVfMjAxOVshaXMubmEoY2l0eV9kYXRhJG1lYW5fcG9wdWxhdGlvbl9leHBvc3VyZV90b19wbTJfNV8yMDE5KV0pCmBgYAoKYGBge3J9CmZldGNoX2NpdHlfZGF0YV9mb3IgPC0gZnVuY3Rpb24ocG9vbF9uYW1lLCBpbmNsdWRlX2NpdHlfbmFtZSA9IEYpIHsKICByZXN1bHRzX2ZpbGVuYW1lIDwtIHBhc3RlKHBhc3RlKCdwZXJjZW50YWdlX29mX3JlZ2lvbmFsX3JpY2huZXNzX19vdXRwdXRfJywgcG9vbF9uYW1lLCAnY2l0eScsICdyaWNobmVzcycsICdpbnRlcmNlcHQnLCBzZXAgPSAiXyIpLCAiY3N2Iiwgc2VwID0gIi4iKQogIHJlc3VsdHMgPC0gcmVhZF9jc3YocmVzdWx0c19maWxlbmFtZSkKICAKICBqb2luZWQgPC0gbGVmdF9qb2luKGNpdHlfZGF0YSwgcmVzdWx0cykKICBqb2luZWQkYWJzX2NpdHlfY2VudHJlX2xhdGl0dWRlID0gYWJzKGpvaW5lZCRjaXR5X2NlbnRyZV9sYXRpdHVkZSkKICAKICByZXF1aXJlZF9jb2x1bW5zIDwtIGMoInBvcHVsYXRpb25fZ3Jvd3RoIiwgInJhaW5mYWxsX21vbnRobHlfbWluIiwgInJhaW5mYWxsX2FubnVhbF9hdmVyYWdlIiwgInJhaW5mYWxsX21vbnRobHlfbWF4IiwgInRlbXBlcmF0dXJlX2FubnVhbF9hdmVyYWdlIiwgInRlbXBlcmF0dXJlX21vbnRobHlfbWluIiwgInRlbXBlcmF0dXJlX21vbnRobHlfbWF4IiwgImhhcHBpbmVzc19uZWdhdGl2ZV9lZmZlY3QiLCAiaGFwcGluZXNzX3Bvc2l0aXZlX2VmZmVjdCIsICJoYXBwaW5lc3NfZnV0dXJlX2xpZmUiLCAibnVtYmVyX29mX2Jpb21lcyIsICJyZWFsbSIsICJiaW9tZV9uYW1lIiwgInJlZ2lvbl8yMGttX2luY2x1ZGVzX2VzdHVhcnkiLCAicmVnaW9uXzUwa21faW5jbHVkZXNfZXN0dWFyeSIsICJyZWdpb25fMTAwa21faW5jbHVkZXNfZXN0dWFyeSIsICJjaXR5X2luY2x1ZGVzX2VzdHVhcnkiLCAicmVnaW9uXzIwa21fYXZlcmFnZV9wb3BfZGVuc2l0eSIsICJyZWdpb25fNTBrbV9hdmVyYWdlX3BvcF9kZW5zaXR5IiwgInJlZ2lvbl8xMDBrbV9hdmVyYWdlX3BvcF9kZW5zaXR5IiwgImNpdHlfbWF4X3BvcF9kZW5zaXR5IiwgImNpdHlfYXZlcmFnZV9wb3BfZGVuc2l0eSIsICJtZWFuX3BvcHVsYXRpb25fZXhwb3N1cmVfdG9fcG0yXzVfMjAxOSIsICJyZWdpb25fMjBrbV9jdWx0aXZhdGVkIiwgInJlZ2lvbl8yMGttX3VyYmFuIiwgInJlZ2lvbl81MGttX2N1bHRpdmF0ZWQiLCAicmVnaW9uXzUwa21fdXJiYW4iLCAicmVnaW9uXzEwMGttX2N1bHRpdmF0ZWQiLCAicmVnaW9uXzEwMGttX3VyYmFuIiwgInJlZ2lvbl8yMGttX2VsZXZhdGlvbl9kZWx0YSIsICJyZWdpb25fMjBrbV9tZWFuX2VsZXZhdGlvbiIsICJyZWdpb25fNTBrbV9lbGV2YXRpb25fZGVsdGEiLCAicmVnaW9uXzUwa21fbWVhbl9lbGV2YXRpb24iLCAicmVnaW9uXzEwMGttX2VsZXZhdGlvbl9kZWx0YSIsICJyZWdpb25fMTAwa21fbWVhbl9lbGV2YXRpb24iLCAiY2l0eV9lbGV2YXRpb25fZGVsdGEiLCAiY2l0eV9tZWFuX2VsZXZhdGlvbiIsICJ1cmJhbiIsICJzaHJ1YnMiLCAicGVybWFuZW50X3dhdGVyIiwgIm9wZW5fZm9yZXN0IiwgImhlcmJhY2VvdXNfd2V0bGFuZCIsICJoZXJiYWNlb3VzX3ZlZ2V0YXRpb24iLCAiY3VsdGl2YXRlZCIsICJjbG9zZWRfZm9yZXN0IiwgInNoYXJlX29mX3BvcHVsYXRpb25fd2l0aGluXzQwMG1fb2Zfb3Blbl9zcGFjZSIsICJwZXJjZW50YWdlX3VyYmFuX2FyZWFfYXNfc3RyZWV0cyIsICJwZXJjZW50YWdlX3VyYmFuX2FyZWFfYXNfb3Blbl9wdWJsaWNfc3BhY2VzX2FuZF9zdHJlZXRzIiwgInBlcmNlbnRhZ2VfdXJiYW5fYXJlYV9hc19vcGVuX3B1YmxpY19zcGFjZXMiLCAiY2l0eV9nZHBfcGVyX3BvcHVsYXRpb24iLCAiY2l0eV9uZHZpIiwgImNpdHlfc3NtIiwgImNpdHlfc3VzbSIsICJyZWdpb25fMjBrbV9uZHZpIiwgInJlZ2lvbl8yMGttX3NzbSIsICJyZWdpb25fMjBrbV9zdXNtIiwgInJlZ2lvbl81MGttX25kdmkiLCAicmVnaW9uXzUwa21fc3NtIiwgInJlZ2lvbl81MGttX3N1c20iLCAicmVnaW9uXzEwMGttX25kdmkiLCAicmVnaW9uXzEwMGttX3NzbSIsICJyZWdpb25fMTAwa21fc3VzbSIsICJjaXR5X3BlcmNlbnRhZ2VfcHJvdGVjdGVkIiwgInJlZ2lvbl8yMGttX3BlcmNlbnRhZ2VfcHJvdGVjdGVkIiwgInJlZ2lvbl81MGttX3BlcmNlbnRhZ2VfcHJvdGVjdGVkIiwgInJlZ2lvbl8xMDBrbV9wZXJjZW50YWdlX3Byb3RlY3RlZCIsICJjaXR5X2NlbnRyZV9sYXRpdHVkZSIsICJhYnNfY2l0eV9jZW50cmVfbGF0aXR1ZGUiKQogIAogIGlmIChpbmNsdWRlX2NpdHlfbmFtZSkgewogICAgcmVxdWlyZWRfY29sdW1ucyA8LSBhcHBlbmQoYygibmFtZSIpLCByZXF1aXJlZF9jb2x1bW5zKQogIH0KICAKICByZXF1aXJlZF9jb2x1bW5zIDwtIGFwcGVuZChjKCJyZXNwb25zZSIpLCByZXF1aXJlZF9jb2x1bW5zKQogIAogIGpvaW5lZFsscmVxdWlyZWRfY29sdW1uc10KfQpgYGAKCgotLS0tLS0tLS0tLS0tLS0tLS0tLS0tCk1lcmxpbiBSZXNwb25zZQotLS0tLS0tLS0tLS0tLS0tLS0tLS0tCgpgYGB7cn0KbWVybGluX2NpdHlfZGF0YSA8LSBmZXRjaF9jaXR5X2RhdGFfZm9yKCdtZXJsaW4nKQptZXJsaW5fY2l0eV9kYXRhCmBgYAoKYGBge3J9Cm1lcmxpbl9jaXR5X2RhdGFfZml4ZWQgPC0gcmZJbXB1dGUocmVzcG9uc2UgfiAuLCBtZXJsaW5fY2l0eV9kYXRhKQptZXJsaW5fY2l0eV9kYXRhX2ZpeGVkCmBgYAoKYGBge3J9CmdncGxvdChtZXJsaW5fY2l0eV9kYXRhX2ZpeGVkLCBhZXMocmVzcG9uc2UpKSArIGdlb21faGlzdG9ncmFtKGJpbndpZHRoID0gMikKYGBgCgpgYGB7cn0KbmFtZXMobWVybGluX2NpdHlfZGF0YV9maXhlZCkKYGBgCgpgYGB7cn0KbWVybGluX3Jlc3BvbnNlIDwtIG1lcmxpbl9jaXR5X2RhdGFfZml4ZWQkcmVzcG9uc2UKbWVybGluX3ByZWRpY3RvcnMgPC0gbWVybGluX2NpdHlfZGF0YV9maXhlZFssLTFdCm1lcmxpbl9wcmVkaWN0b3JzCmBgYAoKYGBge3J9Cm1lcmxpbl9pbnRlcnAgPC0gVlNVUkYoeCA9IG1lcmxpbl9wcmVkaWN0b3JzLCB5ICA9IG1lcmxpbl9yZXNwb25zZSkKYGBgCgpgYGB7cn0KbmFtZXMobWVybGluX3ByZWRpY3RvcnNbLG1lcmxpbl9pbnRlcnAkdmFyc2VsZWN0LmludGVycF0pCmBgYAoKLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQpCaXJkbGlmZSBSZXNwb25zZQotLS0tLS0tLS0tLS0tLS0tLS0tLS0tCmBgYHtyfQpiaXJkbGlmZV9jaXR5X2RhdGEgPC0gZmV0Y2hfY2l0eV9kYXRhX2ZvcignYmlyZGxpZmUnKQpiaXJkbGlmZV9jaXR5X2RhdGEKYGBgCgoKCmBgYHtyfQpnZ3Bsb3QoYmlyZGxpZmVfY2l0eV9kYXRhLCBhZXMocmVzcG9uc2UpKSArIGdlb21faGlzdG9ncmFtKGJpbndpZHRoID0gMSkKYGBgCgpgYGB7cn0KYmlyZGxpZmVfY2l0eV9kYXRhX2ZpeGVkIDwtIHJmSW1wdXRlKHJlc3BvbnNlIH4gLiwgYmlyZGxpZmVfY2l0eV9kYXRhKQpiaXJkbGlmZV9jaXR5X2RhdGFfZml4ZWQKYGBgCgpgYGB7cn0KbmFtZXMoYmlyZGxpZmVfY2l0eV9kYXRhX2ZpeGVkKQpgYGAKCmBgYHtyfQpiaXJkbGlmZV9yZXNwb25zZSA8LSBiaXJkbGlmZV9jaXR5X2RhdGFfZml4ZWQkcmVzcG9uc2UKYmlyZGxpZmVfcHJlZGljdG9ycyA8LSBiaXJkbGlmZV9jaXR5X2RhdGFfZml4ZWRbLC0xXQpiaXJkbGlmZV9wcmVkaWN0b3JzCmBgYAoKCmBgYHtyfQpiaXJkbGlmZV9pbnRlcnAgPC0gVlNVUkYoeCA9IGJpcmRsaWZlX3ByZWRpY3RvcnMsIHkgID0gYmlyZGxpZmVfcmVzcG9uc2UpCmBgYAoKYGBge3J9Cm5hbWVzKGJpcmRsaWZlX3ByZWRpY3RvcnNbLGJpcmRsaWZlX2ludGVycCR2YXJzZWxlY3QuaW50ZXJwXSkKYGBgCgoKLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tClNvLi4uLgotLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KTWVybGluOiAiYWJzX2NpdHlfY2VudHJlX2xhdGl0dWRlIiAgICAgInJlZ2lvbl81MGttX2VsZXZhdGlvbl9kZWx0YSIgICJiaW9tZV9uYW1lIiAgICAgICAgICAgICAgICAgICAicmVnaW9uXzUwa21fc3NtIiAgICAgICAgICAgICAKIFs1XSAicmVnaW9uXzEwMGttX3NzbSIgICAgICAgICAgICAgInJlZ2lvbl8yMGttX2VsZXZhdGlvbl9kZWx0YSIgICJyZWdpb25fMjBrbV91cmJhbiIgICAgICAgICAgICAicmVnaW9uXzEwMGttX2VsZXZhdGlvbl9kZWx0YSIKIFs5XSAidGVtcGVyYXR1cmVfYW5udWFsX2F2ZXJhZ2UiICAgImNpdHlfbmR2aSIgICAgICAgICAgICAgICAgICAgICJjaXR5X2dkcF9wZXJfcG9wdWxhdGlvbiIgICAgICAic2hydWJzIiAgICAgICAgICAgICAgICAgICAgICAKWzEzXSAicGVybWFuZW50X3dhdGVyIiAgICAgICAgICAgICAgImNpdHlfY2VudHJlX2xhdGl0dWRlIiAgICAgICAgICJyZWdpb25fMjBrbV9jdWx0aXZhdGVkIiAgICAgCkJpcmRsaWZlOiAicG9wdWxhdGlvbl9ncm93dGgiICJyZWdpb25fNTBrbV9zc20iICAKCi0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tClRyeSBNb2RlbGxpbmcKLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KCgpgYGB7cn0KbWVybGluX2NpdHlfZGF0YV9uYW1lZCA8LSBmZXRjaF9jaXR5X2RhdGFfZm9yKCdtZXJsaW4nLCBUKQpiaXJkbGlmZV9jaXR5X2RhdGFfbmFtZWQgPC0gZmV0Y2hfY2l0eV9kYXRhX2ZvcignYmlyZGxpZmUnLCBUKQpgYGAKCi0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQpVc2UgY3Jvc3MgdmFsaWRhdGlvbiBhbmQgZHJvcHBpbmcgdGVybXMgdG8gZmluZCBiZXN0IG1vZGVsCi0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQoKZnVsbCBtb2RlbDogIHJlc3BvbnNlIH4gcmVnaW9uXzEwMGttX3NzbSArIHJlZ2lvbl8yMGttX2VsZXZhdGlvbl9kZWx0YSArIHJlZ2lvbl8yMGttX3VyYmFuICsgcmVnaW9uXzEwMGttX2VsZXZhdGlvbl9kZWx0YSArIHRlbXBlcmF0dXJlX2FubnVhbF9hdmVyYWdlICsgY2l0eV9uZHZpICsgY2l0eV9nZHBfcGVyX3BvcHVsYXRpb24gKyBzaHJ1YnMgKyBwZXJtYW5lbnRfd2F0ZXIgKyBjaXR5X2NlbnRyZV9sYXRpdHVkZSArIHJlZ2lvbl8yMGttX2N1bHRpdmF0ZWQKCgpNZXJsaW4gZGF0YSBzZXQKLS0tLS0tLS0tLS0tLS0tLQpgYGB7cn0KbWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwgPC0gbWVybGluX2NpdHlfZGF0YV9maXhlZFttZXJsaW5fY2l0eV9kYXRhX2ZpeGVkJGJpb21lX25hbWUgIT0gJ0JvcmVhbCBGb3Jlc3RzL1RhaWdhJyxdCmBgYAoKYGBge3J9CmN2LmdsbShtZXJsaW5fY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCwgZ2xtKGRhdGEgPSBtZXJsaW5fY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCwgZm9ybXVsYSA9IHJlc3BvbnNlIH4gcmVnaW9uXzEwMGttX3NzbSArIHJlZ2lvbl8yMGttX2VsZXZhdGlvbl9kZWx0YSArIHJlZ2lvbl8yMGttX3VyYmFuICsgcmVnaW9uXzEwMGttX2VsZXZhdGlvbl9kZWx0YSArIHRlbXBlcmF0dXJlX2FubnVhbF9hdmVyYWdlICsgY2l0eV9uZHZpICsgY2l0eV9nZHBfcGVyX3BvcHVsYXRpb24gKyBzaHJ1YnMgKyBwZXJtYW5lbnRfd2F0ZXIgKyBjaXR5X2NlbnRyZV9sYXRpdHVkZSArIHJlZ2lvbl8yMGttX2N1bHRpdmF0ZWQpKSRkZWx0YVsxXQpgYGAKCi0tIENWRSAxOS43Mjg0MQotLSBDYW4gd2UgZHJvcCBvbmU/CgpgYGB7cn0KY3YuZ2xtKG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsLCBnbG0oZm9ybXVsYSA9IHJlc3BvbnNlIH4gcmVnaW9uXzIwa21fZWxldmF0aW9uX2RlbHRhICsgcmVnaW9uXzIwa21fdXJiYW4gKyByZWdpb25fMTAwa21fZWxldmF0aW9uX2RlbHRhICsgdGVtcGVyYXR1cmVfYW5udWFsX2F2ZXJhZ2UgKyBjaXR5X25kdmkgKyBjaXR5X2dkcF9wZXJfcG9wdWxhdGlvbiArIHNocnVicyArIHBlcm1hbmVudF93YXRlciArIGNpdHlfY2VudHJlX2xhdGl0dWRlICsgcmVnaW9uXzIwa21fY3VsdGl2YXRlZCwgZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsKSkkZGVsdGFbMV0KYGBgCgpgYGB7cn0KY3YuZ2xtKG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsLCBnbG0oZm9ybXVsYSA9IHJlc3BvbnNlIH4gcmVnaW9uXzEwMGttX3NzbSArIHJlZ2lvbl8yMGttX3VyYmFuICsgcmVnaW9uXzEwMGttX2VsZXZhdGlvbl9kZWx0YSArIHRlbXBlcmF0dXJlX2FubnVhbF9hdmVyYWdlICsgY2l0eV9uZHZpICsgY2l0eV9nZHBfcGVyX3BvcHVsYXRpb24gKyBzaHJ1YnMgKyBwZXJtYW5lbnRfd2F0ZXIgKyBjaXR5X2NlbnRyZV9sYXRpdHVkZSArIHJlZ2lvbl8yMGttX2N1bHRpdmF0ZWQsIGRhdGEgPSBtZXJsaW5fY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCkpJGRlbHRhWzFdCmBgYAoKYGBge3J9CmN2LmdsbShtZXJsaW5fY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCwgZ2xtKGZvcm11bGEgPSByZXNwb25zZSB+IHJlZ2lvbl8xMDBrbV9zc20gKyByZWdpb25fMjBrbV9lbGV2YXRpb25fZGVsdGEgKyByZWdpb25fMTAwa21fZWxldmF0aW9uX2RlbHRhICsgdGVtcGVyYXR1cmVfYW5udWFsX2F2ZXJhZ2UgKyBjaXR5X25kdmkgKyBjaXR5X2dkcF9wZXJfcG9wdWxhdGlvbiArIHNocnVicyArIHBlcm1hbmVudF93YXRlciArIGNpdHlfY2VudHJlX2xhdGl0dWRlICsgcmVnaW9uXzIwa21fY3VsdGl2YXRlZCwgZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsKSkkZGVsdGFbMV0KYGBgCgpgYGB7cn0KY3YuZ2xtKG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsLCBnbG0oZm9ybXVsYSA9IHJlc3BvbnNlIH4gcmVnaW9uXzEwMGttX3NzbSArIHJlZ2lvbl8yMGttX2VsZXZhdGlvbl9kZWx0YSArIHJlZ2lvbl8yMGttX3VyYmFuICsgdGVtcGVyYXR1cmVfYW5udWFsX2F2ZXJhZ2UgKyBjaXR5X25kdmkgKyBjaXR5X2dkcF9wZXJfcG9wdWxhdGlvbiArIHNocnVicyArIHBlcm1hbmVudF93YXRlciArIGNpdHlfY2VudHJlX2xhdGl0dWRlICsgcmVnaW9uXzIwa21fY3VsdGl2YXRlZCwgZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsKSkkZGVsdGFbMV0KYGBgCgpgYGB7cn0KY3YuZ2xtKG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsLCBnbG0oZm9ybXVsYSA9IHJlc3BvbnNlIH4gcmVnaW9uXzEwMGttX3NzbSArIHJlZ2lvbl8yMGttX2VsZXZhdGlvbl9kZWx0YSArIHJlZ2lvbl8yMGttX3VyYmFuICsgcmVnaW9uXzEwMGttX2VsZXZhdGlvbl9kZWx0YSArIGNpdHlfbmR2aSArIGNpdHlfZ2RwX3Blcl9wb3B1bGF0aW9uICsgc2hydWJzICsgcGVybWFuZW50X3dhdGVyICsgY2l0eV9jZW50cmVfbGF0aXR1ZGUgKyByZWdpb25fMjBrbV9jdWx0aXZhdGVkLCBkYXRhID0gbWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwpKSRkZWx0YVsxXQpgYGAKCmBgYHtyfQpjdi5nbG0obWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwsIGdsbShkYXRhID0gbWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwsIGZvcm11bGEgPSByZXNwb25zZSB+IHJlZ2lvbl8xMDBrbV9zc20gKyByZWdpb25fMjBrbV9lbGV2YXRpb25fZGVsdGEgKyByZWdpb25fMjBrbV91cmJhbiArIHJlZ2lvbl8xMDBrbV9lbGV2YXRpb25fZGVsdGEgKyB0ZW1wZXJhdHVyZV9hbm51YWxfYXZlcmFnZSArIGNpdHlfZ2RwX3Blcl9wb3B1bGF0aW9uICsgc2hydWJzICsgcGVybWFuZW50X3dhdGVyICsgY2l0eV9jZW50cmVfbGF0aXR1ZGUgKyByZWdpb25fMjBrbV9jdWx0aXZhdGVkKSkkZGVsdGFbMV0KYGBgCgpgYGB7cn0KY3YuZ2xtKG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsLCBnbG0oZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsLCBmb3JtdWxhID0gcmVzcG9uc2UgfiByZWdpb25fMTAwa21fc3NtICsgcmVnaW9uXzIwa21fZWxldmF0aW9uX2RlbHRhICsgcmVnaW9uXzIwa21fdXJiYW4gKyByZWdpb25fMTAwa21fZWxldmF0aW9uX2RlbHRhICsgdGVtcGVyYXR1cmVfYW5udWFsX2F2ZXJhZ2UgKyBjaXR5X25kdmkgKyBzaHJ1YnMgKyBwZXJtYW5lbnRfd2F0ZXIgKyBjaXR5X2NlbnRyZV9sYXRpdHVkZSArIHJlZ2lvbl8yMGttX2N1bHRpdmF0ZWQpKSRkZWx0YVsxXQpgYGAKCmBgYHtyfQpjdi5nbG0obWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwsIGdsbShkYXRhID0gbWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwsIGZvcm11bGEgPSByZXNwb25zZSB+IHJlZ2lvbl8xMDBrbV9zc20gKyByZWdpb25fMjBrbV9lbGV2YXRpb25fZGVsdGEgKyByZWdpb25fMjBrbV91cmJhbiArIHJlZ2lvbl8xMDBrbV9lbGV2YXRpb25fZGVsdGEgKyB0ZW1wZXJhdHVyZV9hbm51YWxfYXZlcmFnZSArIGNpdHlfbmR2aSArIGNpdHlfZ2RwX3Blcl9wb3B1bGF0aW9uICsgcGVybWFuZW50X3dhdGVyICsgY2l0eV9jZW50cmVfbGF0aXR1ZGUgKyByZWdpb25fMjBrbV9jdWx0aXZhdGVkKSkkZGVsdGFbMV0KYGBgCgpgYGB7cn0KY3YuZ2xtKG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsLCBnbG0oZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsLCBmb3JtdWxhID0gcmVzcG9uc2UgfiByZWdpb25fMTAwa21fc3NtICsgcmVnaW9uXzIwa21fZWxldmF0aW9uX2RlbHRhICsgcmVnaW9uXzIwa21fdXJiYW4gKyByZWdpb25fMTAwa21fZWxldmF0aW9uX2RlbHRhICsgdGVtcGVyYXR1cmVfYW5udWFsX2F2ZXJhZ2UgKyBjaXR5X25kdmkgKyBjaXR5X2dkcF9wZXJfcG9wdWxhdGlvbiArIHNocnVicyArIGNpdHlfY2VudHJlX2xhdGl0dWRlICsgcmVnaW9uXzIwa21fY3VsdGl2YXRlZCkpJGRlbHRhWzFdCmBgYAoKYGBge3J9CmN2LmdsbShtZXJsaW5fY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCwgZ2xtKGRhdGEgPSBtZXJsaW5fY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCwgZm9ybXVsYSA9IHJlc3BvbnNlIH4gcmVnaW9uXzEwMGttX3NzbSArIHJlZ2lvbl8yMGttX2VsZXZhdGlvbl9kZWx0YSArIHJlZ2lvbl8yMGttX3VyYmFuICsgcmVnaW9uXzEwMGttX2VsZXZhdGlvbl9kZWx0YSArIHRlbXBlcmF0dXJlX2FubnVhbF9hdmVyYWdlICsgY2l0eV9uZHZpICsgY2l0eV9nZHBfcGVyX3BvcHVsYXRpb24gKyBzaHJ1YnMgKyBwZXJtYW5lbnRfd2F0ZXIgKyByZWdpb25fMjBrbV9jdWx0aXZhdGVkKSkkZGVsdGFbMV0KYGBgCgpgYGB7cn0KY3YuZ2xtKG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsLCBnbG0oZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsLCBmb3JtdWxhID0gcmVzcG9uc2UgfiByZWdpb25fMTAwa21fc3NtICsgcmVnaW9uXzIwa21fZWxldmF0aW9uX2RlbHRhICsgcmVnaW9uXzIwa21fdXJiYW4gKyByZWdpb25fMTAwa21fZWxldmF0aW9uX2RlbHRhICsgdGVtcGVyYXR1cmVfYW5udWFsX2F2ZXJhZ2UgKyBjaXR5X25kdmkgKyBjaXR5X2dkcF9wZXJfcG9wdWxhdGlvbiArIHNocnVicyArIHBlcm1hbmVudF93YXRlciArIGNpdHlfY2VudHJlX2xhdGl0dWRlKSkkZGVsdGFbMV0KYGBgCgotLSBkcm9wIGNpdHlfbmR2aSB0byBnaXZlIHNtYWxsZXIgQ1ZFIG9mIDE5LjM1Ci0tIGNhbiB3ZSBkcm9wIGFub3RoZXI/CgpgYGB7cn0KY3YuZ2xtKG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsLCBnbG0oZm9ybXVsYSA9IHJlc3BvbnNlIH4gcmVnaW9uXzIwa21fZWxldmF0aW9uX2RlbHRhICsgcmVnaW9uXzIwa21fdXJiYW4gKyByZWdpb25fMTAwa21fZWxldmF0aW9uX2RlbHRhICsgdGVtcGVyYXR1cmVfYW5udWFsX2F2ZXJhZ2UgKyBjaXR5X2dkcF9wZXJfcG9wdWxhdGlvbiArIHNocnVicyArIHBlcm1hbmVudF93YXRlciArIGNpdHlfY2VudHJlX2xhdGl0dWRlICsgcmVnaW9uXzIwa21fY3VsdGl2YXRlZCwgZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsKSkkZGVsdGFbMV0KYGBgCmBgYHtyfQpjdi5nbG0obWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwsIGdsbShmb3JtdWxhID0gcmVzcG9uc2UgfiByZWdpb25fMTAwa21fc3NtICsgcmVnaW9uXzIwa21fdXJiYW4gKyByZWdpb25fMTAwa21fZWxldmF0aW9uX2RlbHRhICsgdGVtcGVyYXR1cmVfYW5udWFsX2F2ZXJhZ2UgKyBjaXR5X2dkcF9wZXJfcG9wdWxhdGlvbiArIHNocnVicyArIHBlcm1hbmVudF93YXRlciArIGNpdHlfY2VudHJlX2xhdGl0dWRlICsgcmVnaW9uXzIwa21fY3VsdGl2YXRlZCwgZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsKSkkZGVsdGFbMV0KYGBgCgpgYGB7cn0KY3YuZ2xtKG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsLCBnbG0oZm9ybXVsYSA9IHJlc3BvbnNlIH4gcmVnaW9uXzEwMGttX3NzbSArIHJlZ2lvbl8yMGttX2VsZXZhdGlvbl9kZWx0YSArIHJlZ2lvbl8xMDBrbV9lbGV2YXRpb25fZGVsdGEgKyB0ZW1wZXJhdHVyZV9hbm51YWxfYXZlcmFnZSArIGNpdHlfZ2RwX3Blcl9wb3B1bGF0aW9uICsgc2hydWJzICsgcGVybWFuZW50X3dhdGVyICsgY2l0eV9jZW50cmVfbGF0aXR1ZGUgKyByZWdpb25fMjBrbV9jdWx0aXZhdGVkLCBkYXRhID0gbWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwpKSRkZWx0YVsxXQpgYGAKCmBgYHtyfQpjdi5nbG0obWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwsIGdsbShmb3JtdWxhID0gcmVzcG9uc2UgfiByZWdpb25fMTAwa21fc3NtICsgcmVnaW9uXzIwa21fZWxldmF0aW9uX2RlbHRhICsgcmVnaW9uXzIwa21fdXJiYW4gKyByZWdpb25fMTAwa21fZWxldmF0aW9uX2RlbHRhICsgY2l0eV9nZHBfcGVyX3BvcHVsYXRpb24gKyBzaHJ1YnMgKyBwZXJtYW5lbnRfd2F0ZXIgKyBjaXR5X2NlbnRyZV9sYXRpdHVkZSArIHJlZ2lvbl8yMGttX2N1bHRpdmF0ZWQsIGRhdGEgPSBtZXJsaW5fY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCkpJGRlbHRhWzFdCmBgYAoKYGBge3J9CmN2LmdsbShtZXJsaW5fY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCwgZ2xtKGZvcm11bGEgPSByZXNwb25zZSB+IHJlZ2lvbl8xMDBrbV9zc20gKyByZWdpb25fMjBrbV9lbGV2YXRpb25fZGVsdGEgKyByZWdpb25fMjBrbV91cmJhbiArIHJlZ2lvbl8xMDBrbV9lbGV2YXRpb25fZGVsdGEgKyB0ZW1wZXJhdHVyZV9hbm51YWxfYXZlcmFnZSArIHNocnVicyArIHBlcm1hbmVudF93YXRlciArIGNpdHlfY2VudHJlX2xhdGl0dWRlICsgcmVnaW9uXzIwa21fY3VsdGl2YXRlZCwgZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsKSkkZGVsdGFbMV0KYGBgCmBgYHtyfQpjdi5nbG0obWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwsIGdsbShmb3JtdWxhID0gcmVzcG9uc2UgfiByZWdpb25fMTAwa21fc3NtICsgcmVnaW9uXzIwa21fZWxldmF0aW9uX2RlbHRhICsgcmVnaW9uXzIwa21fdXJiYW4gKyByZWdpb25fMTAwa21fZWxldmF0aW9uX2RlbHRhICsgdGVtcGVyYXR1cmVfYW5udWFsX2F2ZXJhZ2UgKyBjaXR5X2dkcF9wZXJfcG9wdWxhdGlvbiArIHBlcm1hbmVudF93YXRlciArIGNpdHlfY2VudHJlX2xhdGl0dWRlICsgcmVnaW9uXzIwa21fY3VsdGl2YXRlZCwgZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsKSkkZGVsdGFbMV0KYGBgCmBgYHtyfQpjdi5nbG0obWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwsIGdsbShmb3JtdWxhID0gcmVzcG9uc2UgfiByZWdpb25fMTAwa21fc3NtICsgcmVnaW9uXzIwa21fZWxldmF0aW9uX2RlbHRhICsgcmVnaW9uXzIwa21fdXJiYW4gKyByZWdpb25fMTAwa21fZWxldmF0aW9uX2RlbHRhICsgdGVtcGVyYXR1cmVfYW5udWFsX2F2ZXJhZ2UgKyBjaXR5X2dkcF9wZXJfcG9wdWxhdGlvbiArIHNocnVicyArIGNpdHlfY2VudHJlX2xhdGl0dWRlICsgcmVnaW9uXzIwa21fY3VsdGl2YXRlZCwgZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsKSkkZGVsdGFbMV0KYGBgCmBgYHtyfQpjdi5nbG0obWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwsIGdsbShmb3JtdWxhID0gcmVzcG9uc2UgfiByZWdpb25fMTAwa21fc3NtICsgcmVnaW9uXzIwa21fZWxldmF0aW9uX2RlbHRhICsgcmVnaW9uXzIwa21fdXJiYW4gKyByZWdpb25fMTAwa21fZWxldmF0aW9uX2RlbHRhICsgdGVtcGVyYXR1cmVfYW5udWFsX2F2ZXJhZ2UgKyBjaXR5X2dkcF9wZXJfcG9wdWxhdGlvbiArIHNocnVicyArIHBlcm1hbmVudF93YXRlciArIHJlZ2lvbl8yMGttX2N1bHRpdmF0ZWQsIGRhdGEgPSBtZXJsaW5fY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCkpJGRlbHRhWzFdCmBgYApgYGB7cn0KY3YuZ2xtKG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsLCBnbG0oZm9ybXVsYSA9IHJlc3BvbnNlIH4gcmVnaW9uXzEwMGttX3NzbSArIHJlZ2lvbl8yMGttX2VsZXZhdGlvbl9kZWx0YSArIHJlZ2lvbl8yMGttX3VyYmFuICsgcmVnaW9uXzEwMGttX2VsZXZhdGlvbl9kZWx0YSArIHRlbXBlcmF0dXJlX2FubnVhbF9hdmVyYWdlICsgY2l0eV9nZHBfcGVyX3BvcHVsYXRpb24gKyBzaHJ1YnMgKyBwZXJtYW5lbnRfd2F0ZXIgKyBjaXR5X2NlbnRyZV9sYXRpdHVkZSwgZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsKSkkZGVsdGFbMV0KYGBgCi0tIGRyb3AgY2l0eV9jZW50cmVfbGF0aXR1ZGUgdG8gZ2l2ZSBzbWFsbGVyIENWRSBvZiAxOS4wNgotLSBjYW4gd2UgZHJvcCBhbm90aGVyPwoKYGBge3J9CmN2LmdsbShtZXJsaW5fY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCwgZ2xtKGZvcm11bGEgPSByZXNwb25zZSB+IHJlZ2lvbl8yMGttX2VsZXZhdGlvbl9kZWx0YSArIHJlZ2lvbl8yMGttX3VyYmFuICsgcmVnaW9uXzEwMGttX2VsZXZhdGlvbl9kZWx0YSArIHRlbXBlcmF0dXJlX2FubnVhbF9hdmVyYWdlICsgY2l0eV9nZHBfcGVyX3BvcHVsYXRpb24gKyBzaHJ1YnMgKyBwZXJtYW5lbnRfd2F0ZXIgKyByZWdpb25fMjBrbV9jdWx0aXZhdGVkLCBkYXRhID0gbWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwpKSRkZWx0YVsxXQpgYGAKCmBgYHtyfQpjdi5nbG0obWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwsIGdsbShmb3JtdWxhID0gcmVzcG9uc2UgfiByZWdpb25fMTAwa21fc3NtICsgcmVnaW9uXzIwa21fdXJiYW4gKyByZWdpb25fMTAwa21fZWxldmF0aW9uX2RlbHRhICsgdGVtcGVyYXR1cmVfYW5udWFsX2F2ZXJhZ2UgKyBjaXR5X2dkcF9wZXJfcG9wdWxhdGlvbiArIHNocnVicyArIHBlcm1hbmVudF93YXRlciArIHJlZ2lvbl8yMGttX2N1bHRpdmF0ZWQsIGRhdGEgPSBtZXJsaW5fY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCkpJGRlbHRhWzFdCmBgYApgYGB7cn0KY3YuZ2xtKG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsLCBnbG0oZm9ybXVsYSA9IHJlc3BvbnNlIH4gcmVnaW9uXzEwMGttX3NzbSArIHJlZ2lvbl8yMGttX2VsZXZhdGlvbl9kZWx0YSArIHJlZ2lvbl8xMDBrbV9lbGV2YXRpb25fZGVsdGEgKyB0ZW1wZXJhdHVyZV9hbm51YWxfYXZlcmFnZSArIGNpdHlfZ2RwX3Blcl9wb3B1bGF0aW9uICsgc2hydWJzICsgcGVybWFuZW50X3dhdGVyICsgcmVnaW9uXzIwa21fY3VsdGl2YXRlZCwgZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsKSkkZGVsdGFbMV0KYGBgCmBgYHtyfQpjdi5nbG0obWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwsIGdsbShmb3JtdWxhID0gcmVzcG9uc2UgfiByZWdpb25fMTAwa21fc3NtICsgcmVnaW9uXzIwa21fZWxldmF0aW9uX2RlbHRhICsgcmVnaW9uXzIwa21fdXJiYW4gKyB0ZW1wZXJhdHVyZV9hbm51YWxfYXZlcmFnZSArIGNpdHlfZ2RwX3Blcl9wb3B1bGF0aW9uICsgc2hydWJzICsgcGVybWFuZW50X3dhdGVyICsgcmVnaW9uXzIwa21fY3VsdGl2YXRlZCwgZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsKSkkZGVsdGFbMV0KYGBgCgpgYGB7cn0KY3YuZ2xtKG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsLCBnbG0oZm9ybXVsYSA9IHJlc3BvbnNlIH4gcmVnaW9uXzEwMGttX3NzbSArIHJlZ2lvbl8yMGttX2VsZXZhdGlvbl9kZWx0YSArIHJlZ2lvbl8yMGttX3VyYmFuICsgcmVnaW9uXzEwMGttX2VsZXZhdGlvbl9kZWx0YSArIGNpdHlfZ2RwX3Blcl9wb3B1bGF0aW9uICsgc2hydWJzICsgcGVybWFuZW50X3dhdGVyICsgcmVnaW9uXzIwa21fY3VsdGl2YXRlZCwgZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsKSkkZGVsdGFbMV0KYGBgCmBgYHtyfQpjdi5nbG0obWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwsIGdsbShmb3JtdWxhID0gcmVzcG9uc2UgfiByZWdpb25fMTAwa21fc3NtICsgcmVnaW9uXzIwa21fZWxldmF0aW9uX2RlbHRhICsgcmVnaW9uXzIwa21fdXJiYW4gKyByZWdpb25fMTAwa21fZWxldmF0aW9uX2RlbHRhICsgdGVtcGVyYXR1cmVfYW5udWFsX2F2ZXJhZ2UgKyBzaHJ1YnMgKyBwZXJtYW5lbnRfd2F0ZXIgKyByZWdpb25fMjBrbV9jdWx0aXZhdGVkLCBkYXRhID0gbWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwpKSRkZWx0YVsxXQpgYGAKYGBge3J9CmN2LmdsbShtZXJsaW5fY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCwgZ2xtKGZvcm11bGEgPSByZXNwb25zZSB+IHJlZ2lvbl8xMDBrbV9zc20gKyByZWdpb25fMjBrbV9lbGV2YXRpb25fZGVsdGEgKyByZWdpb25fMjBrbV91cmJhbiArIHJlZ2lvbl8xMDBrbV9lbGV2YXRpb25fZGVsdGEgKyB0ZW1wZXJhdHVyZV9hbm51YWxfYXZlcmFnZSArIGNpdHlfZ2RwX3Blcl9wb3B1bGF0aW9uICsgcGVybWFuZW50X3dhdGVyICsgcmVnaW9uXzIwa21fY3VsdGl2YXRlZCwgZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsKSkkZGVsdGFbMV0KYGBgCmBgYHtyfQpjdi5nbG0obWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwsIGdsbShmb3JtdWxhID0gcmVzcG9uc2UgfiByZWdpb25fMTAwa21fc3NtICsgcmVnaW9uXzIwa21fZWxldmF0aW9uX2RlbHRhICsgcmVnaW9uXzIwa21fdXJiYW4gKyByZWdpb25fMTAwa21fZWxldmF0aW9uX2RlbHRhICsgdGVtcGVyYXR1cmVfYW5udWFsX2F2ZXJhZ2UgKyBjaXR5X2dkcF9wZXJfcG9wdWxhdGlvbiArIHNocnVicyArIHJlZ2lvbl8yMGttX2N1bHRpdmF0ZWQsIGRhdGEgPSBtZXJsaW5fY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCkpJGRlbHRhWzFdCmBgYAoKYGBge3J9CmN2LmdsbShtZXJsaW5fY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCwgZ2xtKGZvcm11bGEgPSByZXNwb25zZSB+IHJlZ2lvbl8xMDBrbV9zc20gKyByZWdpb25fMjBrbV9lbGV2YXRpb25fZGVsdGEgKyByZWdpb25fMjBrbV91cmJhbiArIHJlZ2lvbl8xMDBrbV9lbGV2YXRpb25fZGVsdGEgKyB0ZW1wZXJhdHVyZV9hbm51YWxfYXZlcmFnZSArIGNpdHlfZ2RwX3Blcl9wb3B1bGF0aW9uICsgc2hydWJzICsgcGVybWFuZW50X3dhdGVyLCBkYXRhID0gbWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwpKSRkZWx0YVsxXQpgYGAKCi0tIGRyb3AgcmVnaW9uXzIwa21fZWxldmF0aW9uX2RlbHRhIHRvIGdpdmUgc21hbGxlciBDVkUgb2YgMTguNzYKLS0gY2FuIHdlIGRyb3AgYW5vdGhlcj8KCmBgYHtyfQpjdi5nbG0obWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwsIGdsbShmb3JtdWxhID0gcmVzcG9uc2UgfiByZWdpb25fMjBrbV91cmJhbiArIHJlZ2lvbl8xMDBrbV9lbGV2YXRpb25fZGVsdGEgKyB0ZW1wZXJhdHVyZV9hbm51YWxfYXZlcmFnZSArIGNpdHlfZ2RwX3Blcl9wb3B1bGF0aW9uICsgc2hydWJzICsgcGVybWFuZW50X3dhdGVyICsgcmVnaW9uXzIwa21fY3VsdGl2YXRlZCwgZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsKSkkZGVsdGFbMV0KYGBgCmBgYHtyfQpjdi5nbG0obWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwsIGdsbShmb3JtdWxhID0gcmVzcG9uc2UgfiByZWdpb25fMTAwa21fc3NtICsgcmVnaW9uXzEwMGttX2VsZXZhdGlvbl9kZWx0YSArIHRlbXBlcmF0dXJlX2FubnVhbF9hdmVyYWdlICsgY2l0eV9nZHBfcGVyX3BvcHVsYXRpb24gKyBzaHJ1YnMgKyBwZXJtYW5lbnRfd2F0ZXIgKyByZWdpb25fMjBrbV9jdWx0aXZhdGVkLCBkYXRhID0gbWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwpKSRkZWx0YVsxXQpgYGAKYGBge3J9CmN2LmdsbShtZXJsaW5fY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCwgZ2xtKGZvcm11bGEgPSByZXNwb25zZSB+IHJlZ2lvbl8xMDBrbV9zc20gKyByZWdpb25fMjBrbV91cmJhbiArIHRlbXBlcmF0dXJlX2FubnVhbF9hdmVyYWdlICsgY2l0eV9nZHBfcGVyX3BvcHVsYXRpb24gKyBzaHJ1YnMgKyBwZXJtYW5lbnRfd2F0ZXIgKyByZWdpb25fMjBrbV9jdWx0aXZhdGVkLCBkYXRhID0gbWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwpKSRkZWx0YVsxXQpgYGAKYGBge3J9CmN2LmdsbShtZXJsaW5fY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCwgZ2xtKGZvcm11bGEgPSByZXNwb25zZSB+IHJlZ2lvbl8xMDBrbV9zc20gKyByZWdpb25fMjBrbV91cmJhbiArIHJlZ2lvbl8xMDBrbV9lbGV2YXRpb25fZGVsdGEgKyBjaXR5X2dkcF9wZXJfcG9wdWxhdGlvbiArIHNocnVicyArIHBlcm1hbmVudF93YXRlciArIHJlZ2lvbl8yMGttX2N1bHRpdmF0ZWQsIGRhdGEgPSBtZXJsaW5fY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCkpJGRlbHRhWzFdCmBgYApgYGB7cn0KY3YuZ2xtKG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsLCBnbG0oZm9ybXVsYSA9IHJlc3BvbnNlIH4gcmVnaW9uXzEwMGttX3NzbSArIHJlZ2lvbl8yMGttX3VyYmFuICsgcmVnaW9uXzEwMGttX2VsZXZhdGlvbl9kZWx0YSArIHRlbXBlcmF0dXJlX2FubnVhbF9hdmVyYWdlICsgc2hydWJzICsgcGVybWFuZW50X3dhdGVyICsgcmVnaW9uXzIwa21fY3VsdGl2YXRlZCwgZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsKSkkZGVsdGFbMV0KYGBgCmBgYHtyfQpjdi5nbG0obWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwsIGdsbShmb3JtdWxhID0gcmVzcG9uc2UgfiByZWdpb25fMTAwa21fc3NtICsgcmVnaW9uXzIwa21fdXJiYW4gKyByZWdpb25fMTAwa21fZWxldmF0aW9uX2RlbHRhICsgdGVtcGVyYXR1cmVfYW5udWFsX2F2ZXJhZ2UgKyBjaXR5X2dkcF9wZXJfcG9wdWxhdGlvbiArIHBlcm1hbmVudF93YXRlciArIHJlZ2lvbl8yMGttX2N1bHRpdmF0ZWQsIGRhdGEgPSBtZXJsaW5fY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCkpJGRlbHRhWzFdCmBgYApgYGB7cn0KY3YuZ2xtKG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsLCBnbG0oZm9ybXVsYSA9IHJlc3BvbnNlIH4gcmVnaW9uXzEwMGttX3NzbSArIHJlZ2lvbl8yMGttX3VyYmFuICsgcmVnaW9uXzEwMGttX2VsZXZhdGlvbl9kZWx0YSArIHRlbXBlcmF0dXJlX2FubnVhbF9hdmVyYWdlICsgY2l0eV9nZHBfcGVyX3BvcHVsYXRpb24gKyBzaHJ1YnMgKyByZWdpb25fMjBrbV9jdWx0aXZhdGVkLCBkYXRhID0gbWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwpKSRkZWx0YVsxXQpgYGAKCmBgYHtyfQpjdi5nbG0obWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwsIGdsbShmb3JtdWxhID0gcmVzcG9uc2UgfiByZWdpb25fMTAwa21fc3NtICsgcmVnaW9uXzIwa21fdXJiYW4gKyByZWdpb25fMTAwa21fZWxldmF0aW9uX2RlbHRhICsgdGVtcGVyYXR1cmVfYW5udWFsX2F2ZXJhZ2UgKyBjaXR5X2dkcF9wZXJfcG9wdWxhdGlvbiArIHNocnVicyArIHBlcm1hbmVudF93YXRlciwgZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsKSkkZGVsdGFbMV0KYGBgCgotLSBkcm9wIHJlZ2lvbl8yMGttX2N1bHRpdmF0ZWQgdG8gZ2l2ZSBzbWFsbGVyIENWRSBvZiAxOC41NAotLSBjYW4gd2UgZHJvcCBhbm90aGVyPwpgYGB7cn0KY3YuZ2xtKG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsLCBnbG0oZm9ybXVsYSA9IHJlc3BvbnNlIH4gcmVnaW9uXzIwa21fdXJiYW4gKyByZWdpb25fMTAwa21fZWxldmF0aW9uX2RlbHRhICsgdGVtcGVyYXR1cmVfYW5udWFsX2F2ZXJhZ2UgKyBjaXR5X2dkcF9wZXJfcG9wdWxhdGlvbiArIHNocnVicyArIHBlcm1hbmVudF93YXRlciwgZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsKSkkZGVsdGFbMV0KYGBgCmBgYHtyfQpjdi5nbG0obWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwsIGdsbShmb3JtdWxhID0gcmVzcG9uc2UgfiByZWdpb25fMTAwa21fc3NtICsgcmVnaW9uXzEwMGttX2VsZXZhdGlvbl9kZWx0YSArIHRlbXBlcmF0dXJlX2FubnVhbF9hdmVyYWdlICsgY2l0eV9nZHBfcGVyX3BvcHVsYXRpb24gKyBzaHJ1YnMgKyBwZXJtYW5lbnRfd2F0ZXIsIGRhdGEgPSBtZXJsaW5fY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCkpJGRlbHRhWzFdCmBgYApgYGB7cn0KY3YuZ2xtKG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsLCBnbG0oZm9ybXVsYSA9IHJlc3BvbnNlIH4gcmVnaW9uXzEwMGttX3NzbSArIHJlZ2lvbl8yMGttX3VyYmFuICsgdGVtcGVyYXR1cmVfYW5udWFsX2F2ZXJhZ2UgKyBjaXR5X2dkcF9wZXJfcG9wdWxhdGlvbiArIHNocnVicyArIHBlcm1hbmVudF93YXRlciwgZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsKSkkZGVsdGFbMV0KYGBgCgpgYGB7cn0KY3YuZ2xtKG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsLCBnbG0oZm9ybXVsYSA9IHJlc3BvbnNlIH4gcmVnaW9uXzEwMGttX3NzbSArIHJlZ2lvbl8yMGttX3VyYmFuICsgcmVnaW9uXzEwMGttX2VsZXZhdGlvbl9kZWx0YSArIGNpdHlfZ2RwX3Blcl9wb3B1bGF0aW9uICsgc2hydWJzICsgcGVybWFuZW50X3dhdGVyLCBkYXRhID0gbWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwpKSRkZWx0YVsxXQpgYGAKYGBge3J9CmN2LmdsbShtZXJsaW5fY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCwgZ2xtKGZvcm11bGEgPSByZXNwb25zZSB+IHJlZ2lvbl8xMDBrbV9zc20gKyByZWdpb25fMjBrbV91cmJhbiArIHJlZ2lvbl8xMDBrbV9lbGV2YXRpb25fZGVsdGEgKyB0ZW1wZXJhdHVyZV9hbm51YWxfYXZlcmFnZSArIHNocnVicyArIHBlcm1hbmVudF93YXRlciwgZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsKSkkZGVsdGFbMV0KYGBgCmBgYHtyfQpjdi5nbG0obWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwsIGdsbShmb3JtdWxhID0gcmVzcG9uc2UgfiByZWdpb25fMTAwa21fc3NtICsgcmVnaW9uXzIwa21fdXJiYW4gKyByZWdpb25fMTAwa21fZWxldmF0aW9uX2RlbHRhICsgdGVtcGVyYXR1cmVfYW5udWFsX2F2ZXJhZ2UgKyBjaXR5X2dkcF9wZXJfcG9wdWxhdGlvbiArIHBlcm1hbmVudF93YXRlciwgZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsKSkkZGVsdGFbMV0KYGBgCmBgYHtyfQpjdi5nbG0obWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwsIGdsbShmb3JtdWxhID0gcmVzcG9uc2UgfiByZWdpb25fMTAwa21fc3NtICsgcmVnaW9uXzIwa21fdXJiYW4gKyByZWdpb25fMTAwa21fZWxldmF0aW9uX2RlbHRhICsgdGVtcGVyYXR1cmVfYW5udWFsX2F2ZXJhZ2UgKyBjaXR5X2dkcF9wZXJfcG9wdWxhdGlvbiArIHNocnVicywgZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsKSkkZGVsdGFbMV0KYGBgCi0tIGRyb3AgcGVybWFuZW50X3dhdGVyIHRvIGdpdmUgc21hbGxlciBDVkUgb2YgMTguMzQKLS0gY2FuIHdlIGRyb3AgYW5vdGhlcj8KCmBgYHtyfQpjdi5nbG0obWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwsIGdsbShmb3JtdWxhID0gcmVzcG9uc2UgfiByZWdpb25fMjBrbV91cmJhbiArIHJlZ2lvbl8xMDBrbV9lbGV2YXRpb25fZGVsdGEgKyB0ZW1wZXJhdHVyZV9hbm51YWxfYXZlcmFnZSArIGNpdHlfZ2RwX3Blcl9wb3B1bGF0aW9uICsgc2hydWJzLCBkYXRhID0gbWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwpKSRkZWx0YVsxXQpgYGAKYGBge3J9CmN2LmdsbShtZXJsaW5fY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCwgZ2xtKGZvcm11bGEgPSByZXNwb25zZSB+IHJlZ2lvbl8xMDBrbV9zc20gKyByZWdpb25fMTAwa21fZWxldmF0aW9uX2RlbHRhICsgdGVtcGVyYXR1cmVfYW5udWFsX2F2ZXJhZ2UgKyBjaXR5X2dkcF9wZXJfcG9wdWxhdGlvbiArIHNocnVicywgZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsKSkkZGVsdGFbMV0KYGBgCmBgYHtyfQpjdi5nbG0obWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwsIGdsbShmb3JtdWxhID0gcmVzcG9uc2UgfiByZWdpb25fMTAwa21fc3NtICsgcmVnaW9uXzIwa21fdXJiYW4gKyB0ZW1wZXJhdHVyZV9hbm51YWxfYXZlcmFnZSArIGNpdHlfZ2RwX3Blcl9wb3B1bGF0aW9uICsgc2hydWJzLCBkYXRhID0gbWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwpKSRkZWx0YVsxXQpgYGAKYGBge3J9CmN2LmdsbShtZXJsaW5fY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCwgZ2xtKGZvcm11bGEgPSByZXNwb25zZSB+IHJlZ2lvbl8xMDBrbV9zc20gKyByZWdpb25fMjBrbV91cmJhbiArIHJlZ2lvbl8xMDBrbV9lbGV2YXRpb25fZGVsdGEgKyBjaXR5X2dkcF9wZXJfcG9wdWxhdGlvbiArIHNocnVicywgZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsKSkkZGVsdGFbMV0KYGBgCmBgYHtyfQpjdi5nbG0obWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwsIGdsbShmb3JtdWxhID0gcmVzcG9uc2UgfiByZWdpb25fMTAwa21fc3NtICsgcmVnaW9uXzIwa21fdXJiYW4gKyByZWdpb25fMTAwa21fZWxldmF0aW9uX2RlbHRhICsgdGVtcGVyYXR1cmVfYW5udWFsX2F2ZXJhZ2UgKyBzaHJ1YnMsIGRhdGEgPSBtZXJsaW5fY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCkpJGRlbHRhWzFdCmBgYApgYGB7cn0KY3YuZ2xtKG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsLCBnbG0oZm9ybXVsYSA9IHJlc3BvbnNlIH4gcmVnaW9uXzEwMGttX3NzbSArIHJlZ2lvbl8yMGttX3VyYmFuICsgcmVnaW9uXzEwMGttX2VsZXZhdGlvbl9kZWx0YSArIHRlbXBlcmF0dXJlX2FubnVhbF9hdmVyYWdlICsgY2l0eV9nZHBfcGVyX3BvcHVsYXRpb24sIGRhdGEgPSBtZXJsaW5fY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCkpJGRlbHRhWzFdCmBgYAoKLS0gZHJvcCB0ZW1wZXJhdHVyZV9hbm51YWxfYXZlcmFnZSB0byBnaXZlIHNtYWxsZXIgQ1ZFIG9mIDE4LjE0Ci0tIGNhbiB3ZSBkcm9wIGFub3RoZXI/CmBgYHtyfQpjdi5nbG0obWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwsIGdsbShmb3JtdWxhID0gcmVzcG9uc2UgfiByZWdpb25fMjBrbV91cmJhbiArIHJlZ2lvbl8xMDBrbV9lbGV2YXRpb25fZGVsdGEgKyBjaXR5X2dkcF9wZXJfcG9wdWxhdGlvbiArIHNocnVicywgZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsKSkkZGVsdGFbMV0KYGBgCgpgYGB7cn0KY3YuZ2xtKG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsLCBnbG0oZm9ybXVsYSA9IHJlc3BvbnNlIH4gcmVnaW9uXzEwMGttX3NzbSArIHJlZ2lvbl8xMDBrbV9lbGV2YXRpb25fZGVsdGEgKyBjaXR5X2dkcF9wZXJfcG9wdWxhdGlvbiArIHNocnVicywgZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsKSkkZGVsdGFbMV0KYGBgCmBgYHtyfQpjdi5nbG0obWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwsIGdsbShmb3JtdWxhID0gcmVzcG9uc2UgfiByZWdpb25fMTAwa21fc3NtICsgcmVnaW9uXzIwa21fdXJiYW4gKyBjaXR5X2dkcF9wZXJfcG9wdWxhdGlvbiArIHNocnVicywgZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsKSkkZGVsdGFbMV0KYGBgCmBgYHtyfQpjdi5nbG0obWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwsIGdsbShmb3JtdWxhID0gcmVzcG9uc2UgfiByZWdpb25fMTAwa21fc3NtICsgcmVnaW9uXzIwa21fdXJiYW4gKyByZWdpb25fMTAwa21fZWxldmF0aW9uX2RlbHRhICsgc2hydWJzLCBkYXRhID0gbWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwpKSRkZWx0YVsxXQpgYGAKYGBge3J9CmN2LmdsbShtZXJsaW5fY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCwgZ2xtKGZvcm11bGEgPSByZXNwb25zZSB+IHJlZ2lvbl8xMDBrbV9zc20gKyByZWdpb25fMjBrbV91cmJhbiArIHJlZ2lvbl8xMDBrbV9lbGV2YXRpb25fZGVsdGEgKyBjaXR5X2dkcF9wZXJfcG9wdWxhdGlvbiwgZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsKSkkZGVsdGFbMV0KYGBgCi0tIGRyb3AgcmVnaW9uXzIwa21fdXJiYW4gdG8gZ2l2ZSBzbWFsbGVyIENWRSBvZiAxOC4wMwotLSBjYW4gd2UgZHJvcCBhbm90aGVyPwoKYGBge3J9CmN2LmdsbShtZXJsaW5fY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCwgZ2xtKGZvcm11bGEgPSByZXNwb25zZSB+IHJlZ2lvbl8xMDBrbV9lbGV2YXRpb25fZGVsdGEgKyBjaXR5X2dkcF9wZXJfcG9wdWxhdGlvbiArIHNocnVicywgZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsKSkkZGVsdGFbMV0KYGBgCgoKYGBge3J9CmN2LmdsbShtZXJsaW5fY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCwgZ2xtKGZvcm11bGEgPSByZXNwb25zZSB+IHJlZ2lvbl8xMDBrbV9zc20gKyBjaXR5X2dkcF9wZXJfcG9wdWxhdGlvbiArIHNocnVicywgZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsKSkkZGVsdGFbMV0KYGBgCgoKYGBge3J9CmN2LmdsbShtZXJsaW5fY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCwgZ2xtKGZvcm11bGEgPSByZXNwb25zZSB+IHJlZ2lvbl8xMDBrbV9zc20gKyByZWdpb25fMTAwa21fZWxldmF0aW9uX2RlbHRhICsgc2hydWJzLCBkYXRhID0gbWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwpKSRkZWx0YVsxXQpgYGAKCmBgYHtyfQpjdi5nbG0obWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwsIGdsbShmb3JtdWxhID0gcmVzcG9uc2UgfiByZWdpb25fMTAwa21fc3NtICsgcmVnaW9uXzEwMGttX2VsZXZhdGlvbl9kZWx0YSArIGNpdHlfZ2RwX3Blcl9wb3B1bGF0aW9uLCBkYXRhID0gbWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwpKSRkZWx0YVsxXQpgYGAKLS0gZHJvcCBzaHJ1YnMgdG8gZ2l2ZSBzbWFsbGVyIENWRSBvZiAxNy45NQotLSBjYW4gd2UgZHJvcCBhbm90aGVyPwoKYGBge3J9CmN2LmdsbShtZXJsaW5fY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCwgZ2xtKGZvcm11bGEgPSByZXNwb25zZSB+IHJlZ2lvbl8xMDBrbV9lbGV2YXRpb25fZGVsdGEgKyBjaXR5X2dkcF9wZXJfcG9wdWxhdGlvbiwgZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsKSkkZGVsdGFbMV0KYGBgCgpgYGB7cn0KY3YuZ2xtKG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsLCBnbG0oZm9ybXVsYSA9IHJlc3BvbnNlIH4gcmVnaW9uXzEwMGttX3NzbSArIGNpdHlfZ2RwX3Blcl9wb3B1bGF0aW9uLCBkYXRhID0gbWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwpKSRkZWx0YVsxXQpgYGAKCmBgYHtyfQpjdi5nbG0obWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwsIGdsbShmb3JtdWxhID0gcmVzcG9uc2UgfiByZWdpb25fMTAwa21fc3NtICsgcmVnaW9uXzEwMGttX2VsZXZhdGlvbl9kZWx0YSwgZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsKSkkZGVsdGFbMV0KYGBgCgotLSBkcm9wIGNpdHlfZ2RwX3Blcl9wb3B1bGF0aW9uIHRvIGdpdmUgc21hbGxlciBDVkUgb2YgMTcuOTQKLS0gY2FuIHdlIGRyb3AgYW5vdGhlcj8KYGBge3J9CmN2LmdsbShtZXJsaW5fY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCwgZ2xtKGZvcm11bGEgPSByZXNwb25zZSB+IHJlZ2lvbl8xMDBrbV9lbGV2YXRpb25fZGVsdGEsIGRhdGEgPSBtZXJsaW5fY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCkpJGRlbHRhWzFdCmBgYAoKYGBge3J9CmN2LmdsbShtZXJsaW5fY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCwgZ2xtKGZvcm11bGEgPSByZXNwb25zZSB+IHJlZ2lvbl8xMDBrbV9zc20sIGRhdGEgPSBtZXJsaW5fY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCkpJGRlbHRhWzFdCmBgYAotLSBObwoKLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQotLSBiZXN0IG1vZGVsIHdpdGggcmVnaW9uXzEwMGttX3NzbSArIHJlZ2lvbl8xMDBrbV9lbGV2YXRpb25fZGVsdGEgKENWIGVycm9yIDE3LjkxMjE2KQotLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tCgpgYGB7cn0Kc3VtbWFyeShnbG0oZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWQsIGZvcm11bGEgPSByZXNwb25zZSB+IHJlZ2lvbl8xMDBrbV9zc20gKyByZWdpb25fMTAwa21fZWxldmF0aW9uX2RlbHRhKSkKYGBgCgpgYGB7cn0KcmVnX21lcmxpbiA9IGdsbShkYXRhID0gbWVybGluX2NpdHlfZGF0YV9maXhlZCwgZm9ybXVsYSA9IHJlc3BvbnNlIH4gcmVnaW9uXzEwMGttX3NzbSArIHJlZ2lvbl8xMDBrbV9lbGV2YXRpb25fZGVsdGEpCndpdGgoc3VtbWFyeShyZWdfbWVybGluKSwgMSAtIGRldmlhbmNlL251bGwuZGV2aWFuY2UpCmBgYAoKQmlyZGxpZmUgZGF0YSBzZXQKLS0tLS0tLS0tLS0tLS0tLQpgYGB7cn0KYmlyZGxpZmVfY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCA8LSBiaXJkbGlmZV9jaXR5X2RhdGFfZml4ZWRbYmlyZGxpZmVfY2l0eV9kYXRhX2ZpeGVkJGJpb21lX25hbWUgIT0gJ0JvcmVhbCBGb3Jlc3RzL1RhaWdhJyxdCmBgYAoKYGBge3J9CmN2LmdsbShiaXJkbGlmZV9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsLCBnbG0oZm9ybXVsYSA9IHJlc3BvbnNlIH4gcmVnaW9uXzEwMGttX3NzbSArIHJlZ2lvbl81MGttX3NzbSArIHJlZ2lvbl81MGttX2VsZXZhdGlvbl9kZWx0YSArIGJpb21lX25hbWUgKyBwb3B1bGF0aW9uX2dyb3d0aCwgZGF0YSA9IGJpcmRsaWZlX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwpKSRkZWx0YVsxXQpgYGAKCi0tIGNhbiB3ZSBkcm9wIGEgdmFyaWFibGU/CgpgYGB7cn0KY3YuZ2xtKGJpcmRsaWZlX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwsIGdsbShmb3JtdWxhID0gcmVzcG9uc2UgfiByZWdpb25fNTBrbV9zc20gKyByZWdpb25fNTBrbV9lbGV2YXRpb25fZGVsdGEgKyBiaW9tZV9uYW1lICsgcG9wdWxhdGlvbl9ncm93dGgsIGRhdGEgPSBiaXJkbGlmZV9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsKSkkZGVsdGFbMV0KYGBgCgpgYGB7cn0KY3YuZ2xtKGJpcmRsaWZlX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwsIGdsbShmb3JtdWxhID0gcmVzcG9uc2UgfiByZWdpb25fMTAwa21fc3NtICsgcmVnaW9uXzUwa21fZWxldmF0aW9uX2RlbHRhICsgYmlvbWVfbmFtZSArIHBvcHVsYXRpb25fZ3Jvd3RoLCBkYXRhID0gYmlyZGxpZmVfY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCkpJGRlbHRhWzFdCmBgYAoKYGBge3J9CmN2LmdsbShiaXJkbGlmZV9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsLCBnbG0oZm9ybXVsYSA9IHJlc3BvbnNlIH4gcmVnaW9uXzEwMGttX3NzbSArIHJlZ2lvbl81MGttX3NzbSArIGJpb21lX25hbWUgKyBwb3B1bGF0aW9uX2dyb3d0aCwgZGF0YSA9IGJpcmRsaWZlX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwpKSRkZWx0YVsxXQpgYGAKCmBgYHtyfQpjdi5nbG0oYmlyZGxpZmVfY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCwgZ2xtKGZvcm11bGEgPSByZXNwb25zZSB+IHJlZ2lvbl8xMDBrbV9zc20gKyByZWdpb25fNTBrbV9zc20gKyByZWdpb25fNTBrbV9lbGV2YXRpb25fZGVsdGEgKyBwb3B1bGF0aW9uX2dyb3d0aCwgZGF0YSA9IGJpcmRsaWZlX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwpKSRkZWx0YVsxXQpgYGAKCmBgYHtyfQpjdi5nbG0oYmlyZGxpZmVfY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCwgZ2xtKGZvcm11bGEgPSByZXNwb25zZSB+IHJlZ2lvbl8xMDBrbV9zc20gKyByZWdpb25fNTBrbV9zc20gKyByZWdpb25fNTBrbV9lbGV2YXRpb25fZGVsdGEgKyBiaW9tZV9uYW1lLCBkYXRhID0gYmlyZGxpZmVfY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCkpJGRlbHRhWzFdCmBgYAoKLS0gZHJvcCBiaW9tZV9uYW1lIHRvIGdpdmUgQ1ZFIG9mIDYuNTAzNDIxCi0tIGNhbiB3ZSBkcm9wIGFub3RoZXI/CgpgYGB7cn0KY3YuZ2xtKGJpcmRsaWZlX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwsIGdsbShmb3JtdWxhID0gcmVzcG9uc2UgfiByZWdpb25fNTBrbV9zc20gKyByZWdpb25fNTBrbV9lbGV2YXRpb25fZGVsdGEgKyBwb3B1bGF0aW9uX2dyb3d0aCwgZGF0YSA9IGJpcmRsaWZlX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwpKSRkZWx0YVsxXQpgYGAKCmBgYHtyfQpjdi5nbG0oYmlyZGxpZmVfY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCwgZ2xtKGZvcm11bGEgPSByZXNwb25zZSB+IHJlZ2lvbl8xMDBrbV9zc20gKyByZWdpb25fNTBrbV9lbGV2YXRpb25fZGVsdGEgKyBwb3B1bGF0aW9uX2dyb3d0aCwgZGF0YSA9IGJpcmRsaWZlX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwpKSRkZWx0YVsxXQpgYGAKCgpgYGB7cn0KY3YuZ2xtKGJpcmRsaWZlX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwsIGdsbShmb3JtdWxhID0gcmVzcG9uc2UgfiByZWdpb25fMTAwa21fc3NtICsgcmVnaW9uXzUwa21fc3NtICsgcG9wdWxhdGlvbl9ncm93dGgsIGRhdGEgPSBiaXJkbGlmZV9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsKSkkZGVsdGFbMV0KYGBgCgoKYGBge3J9CmN2LmdsbShiaXJkbGlmZV9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsLCBnbG0oZm9ybXVsYSA9IHJlc3BvbnNlIH4gcmVnaW9uXzEwMGttX3NzbSArIHJlZ2lvbl81MGttX3NzbSArIHJlZ2lvbl81MGttX2VsZXZhdGlvbl9kZWx0YSwgZGF0YSA9IGJpcmRsaWZlX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwpKSRkZWx0YVsxXQpgYGAKCi0tIGRyb3AgcmVnaW9uXzEwMGttX3NzbSB0byBnaXZlIENWRSBvZiA2LjQxNzMxMQotLSBjYW4gd2UgZHJvcCBhbm90aGVyPwoKYGBge3J9CmN2LmdsbShiaXJkbGlmZV9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsLCBnbG0oZm9ybXVsYSA9IHJlc3BvbnNlIH4gcmVnaW9uXzUwa21fZWxldmF0aW9uX2RlbHRhICsgcG9wdWxhdGlvbl9ncm93dGgsIGRhdGEgPSBiaXJkbGlmZV9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsKSkkZGVsdGFbMV0KYGBgCgpgYGB7cn0KY3YuZ2xtKGJpcmRsaWZlX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwsIGdsbShmb3JtdWxhID0gcmVzcG9uc2UgfiByZWdpb25fNTBrbV9zc20gKyBwb3B1bGF0aW9uX2dyb3d0aCwgZGF0YSA9IGJpcmRsaWZlX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwpKSRkZWx0YVsxXQpgYGAKCmBgYHtyfQpjdi5nbG0oYmlyZGxpZmVfY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCwgZ2xtKGZvcm11bGEgPSByZXNwb25zZSB+IHJlZ2lvbl81MGttX3NzbSArIHJlZ2lvbl81MGttX2VsZXZhdGlvbl9kZWx0YSwgZGF0YSA9IGJpcmRsaWZlX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwpKSRkZWx0YVsxXQpgYGAKCi0tIGRyb3AgcmVnaW9uXzUwa21fZWxldmF0aW9uX2RlbHRhIHRvIGdpdmUgQ1ZFIG9mIDYuMzQyMDI1Ci0tIGNhbiB3ZSBkcm9wIGFub3RoZXI/CgpgYGB7cn0KY3YuZ2xtKGJpcmRsaWZlX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwsIGdsbShmb3JtdWxhID0gcmVzcG9uc2UgfiBwb3B1bGF0aW9uX2dyb3d0aCwgZGF0YSA9IGJpcmRsaWZlX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwpKSRkZWx0YVsxXQpgYGAKCmBgYHtyfQpjdi5nbG0oYmlyZGxpZmVfY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCwgZ2xtKGZvcm11bGEgPSByZXNwb25zZSB+IHJlZ2lvbl81MGttX3NzbSwgZGF0YSA9IGJpcmRsaWZlX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWwpKSRkZWx0YVsxXQpgYGAKCi0tIGRyb3AgY2l0eV9nZHBfcGVyX3BvcHVsYXRpb24gdG8gZ2l2ZSBDVkUgb2YgNi4yOTEyOTkKLS0gaXMgdGhpcyBiZXR0ZXIgdGhhbiBubyB2YXJpYWJsZT8KCmBgYHtyfQpjdi5nbG0oYmlyZGxpZmVfY2l0eV9kYXRhX2ZpeGVkX25vX2JvcmVhbCwgZ2xtKGZvcm11bGEgPSByZXNwb25zZSB+IDEsIGRhdGEgPSBiaXJkbGlmZV9jaXR5X2RhdGFfZml4ZWRfbm9fYm9yZWFsKSkkZGVsdGFbMV0KYGBgCi0tIHllcywganVzdCEKCi0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KLS0gc28gYmVzdCBtb2RlbCB3aXRoIGJpcmRsaWZlIGlzIHJlZ2lvbl81MGttX3NzbQotLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tCgpgYGB7cn0Kc3VtbWFyeShnbG0oZGF0YSA9IGJpcmRsaWZlX2NpdHlfZGF0YV9maXhlZCwgZm9ybXVsYSA9IHJlc3BvbnNlIH4gcmVnaW9uXzUwa21fc3NtKSkKYGBgCgpgYGB7cn0KcmVnX2JpcmRsaWZlID0gZ2xtKGRhdGEgPSBiaXJkbGlmZV9jaXR5X2RhdGFfZml4ZWQsIGZvcm11bGEgPSByZXNwb25zZSB+IHJlZ2lvbl81MGttX3NzbSkKd2l0aChzdW1tYXJ5KHJlZ19iaXJkbGlmZSksIDEgLSBkZXZpYW5jZS9udWxsLmRldmlhbmNlKQpgYGAKCgpgYGB7cn0KZ2dwbG90KGJpcmRsaWZlX2NpdHlfZGF0YV9uYW1lZCwgYWVzKHggPSByZWdpb25fNTBrbV9zc20sIHkgPSByZXNwb25zZSkpICsgZ2VvbV9wb2ludCgpICsgZ2VvbV9zbW9vdGgobWV0aG9kID0gImdsbSIpCmBgYAoKCgoKLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tCkNoZWNrIGJpcmRsaWZlIG1vZGVsIGZpdAotLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KCmBgYHtyfQpiaXJkbGlmZS5maXQgPC0gZ2xtKGRhdGEgPSBiaXJkbGlmZV9jaXR5X2RhdGFfZml4ZWQsIGZvcm11bGEgPSByZXNwb25zZSB+IHJlZ2lvbl81MGttX3NzbSkKc3VtbWFyeShiaXJkbGlmZS5maXQpCndpdGgoc3VtbWFyeShiaXJkbGlmZS5maXQpLCAxIC0gZGV2aWFuY2UvbnVsbC5kZXZpYW5jZSkKcGxvdChiaXJkbGlmZS5maXQpCmBgYAoKYGBge3J9CmJpcmRsaWZlX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWxbYygxNiwgNTMsIDcyKSwgYygicmVnaW9uXzUwa21fc3NtIildCmBgYAoKYGBge3J9CmNpdHlfZGF0YVtjKDE2LCA1MywgNzIpLCBjKCJuYW1lIiwgInJlZ2lvbl81MGttX3NzbSIpXQpgYGAKCmBgYHtyfQpkYXQgPC0gcHJlZGljdChnbG0oZm9ybXVsYSA9IHJlc3BvbnNlIH4gcmVnaW9uXzUwa21fc3NtLCBkYXRhID0gYmlyZGxpZmVfY2l0eV9kYXRhX25hbWVkKSwgc2UuZml0PVQpCgpvdXRzaWRlX3NlIDwtIGJpcmRsaWZlX2NpdHlfZGF0YV9uYW1lZFtiaXJkbGlmZV9jaXR5X2RhdGFfbmFtZWQkcmVzcG9uc2UgPCBkYXQkZml0IC0gMTUqIGRhdCRzZS5maXQgfCBiaXJkbGlmZV9jaXR5X2RhdGFfbmFtZWQkcmVzcG9uc2UgPiBkYXQkZml0ICsgMTUgKiBkYXQkc2UuZml0LF0KCmdncGxvdChiaXJkbGlmZV9jaXR5X2RhdGFfbmFtZWQsIGFlcyh4ID0gcmVnaW9uXzUwa21fc3NtLCB5ID0gcmVzcG9uc2UpKSArIAogIGdlb21fcG9pbnQoc2l6ZT0xKSArIAogIGdlb21fc21vb3RoKG1ldGhvZCA9ICJnbG0iKSArCiAgZ2VvbV90ZXh0KGFlcyhsYWJlbCA9IG5hbWUpLCBkYXRhID0gb3V0c2lkZV9zZSwgc2l6ZSA9IDMsIHBvc2l0aW9uID0gImRvZGdlIiwgdmp1c3QgPSAiaW53YXJkIiwgaGp1c3QgPSAiaW53YXJkIiwgY29sb3IgPSAicmVkIiwgYW5nbGU9LTE1KSArCiAgZ2VvbV9wb2ludChkYXRhID0gb3V0c2lkZV9zZSwgY29sb3IgPSAicmVkIikgKwogIHRoZW1lX2J3KCkgKwogIHlsYWIoIkNpdHkgUmFuZG9tIEVmZmVjdCBJbnRlcmNlcHQiKSArIHhsYWIoIlJlZ2lvbmFsICg1MGttKSBTU00iKSArIGxhYnModGl0bGUgPSAiQmlyZGxpZmUiKQoKZ2dzYXZlKCJjaXR5X2VmZmVjdF9yaWNobmVzc19fb3V0cHV0X19iaXJkbGlmZS5qcGciKQpgYGAKCkhvdyBtdWNoIHZhcmlhdGlvbiBoYXZlIHdlIGV4cGxhaW5lZD8KLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tCgpgYGB7cn0KYmlyZGxpZmVfY2l0eV9kYXRhX25hbWVkJHJlc2lkdWFscyA8LSByZXNpZChiaXJkbGlmZS5maXQpCmdncGxvdChiaXJkbGlmZV9jaXR5X2RhdGFfbmFtZWQsIGFlcyh5ID0gcmVzcG9uc2UsIHggPSByZXNpZHVhbHMpKSArIAogIGdlb21fcG9pbnQoKSArCiAgZ2VvbV9zbW9vdGgobWV0aG9kID0gImxtIiwgc2UgPSBGKSArCiAgbGFicyh0aXRsZSA9ICJCaXJkbGlmZSIsIHN1YnRpdGxlID0gcGFzdGUoIkNvcnJlbGF0aW9uIiwgY29yKGJpcmRsaWZlX2NpdHlfZGF0YV9uYW1lZCRyZXNpZHVhbHMsIGJpcmRsaWZlX2NpdHlfZGF0YV9uYW1lZCRyZXNwb25zZSkpKSArCiAgdGhlbWVfYncoKQpgYGAKCi0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQpDaGVjayBNZXJsaW4gbW9kZWwgZml0Ci0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQoKYGBge3J9Cm1lcmxpbi5maXQgPC0gZ2xtKGRhdGEgPSBtZXJsaW5fY2l0eV9kYXRhX25hbWVkLCBmb3JtdWxhID0gcmVzcG9uc2UgfiByZWdpb25fMTAwa21fc3NtICsgcmVnaW9uXzUwa21fZWxldmF0aW9uX2RlbHRhKQpzdW1tYXJ5KG1lcmxpbi5maXQpCndpdGgoc3VtbWFyeShtZXJsaW4uZml0KSwgMSAtIGRldmlhbmNlL251bGwuZGV2aWFuY2UpCnBsb3QobWVybGluLmZpdCkKYGBgCgpgYGB7cn0KbWVybGluX2NpdHlfZGF0YV9maXhlZF9ub19ib3JlYWxbYygyNCwgNDIsIDcyKSwgYygicmVnaW9uXzEwMGttX3NzbSIsICJyZWdpb25fNTBrbV9lbGV2YXRpb25fZGVsdGEiKV0KYGBgCgpgYGB7cn0KY2l0eV9kYXRhW2MoMjQsIDQyLCA3MiksIGMoIm5hbWUiLCAicmVnaW9uXzEwMGttX3NzbSIsICJyZWdpb25fNTBrbV9lbGV2YXRpb25fZGVsdGEiKV0KYGBgCgpgYGB7cn0KZ2dwbG90KG1lcmxpbl9jaXR5X2RhdGFfbmFtZWQsIGFlcyh4ID0gcmVnaW9uXzEwMGttX3NzbSwgeSA9IHJlc3BvbnNlKSkgKyAKICBnZW9tX3BvaW50KGFlcyhzaXplID0gcmVnaW9uXzUwa21fZWxldmF0aW9uX2RlbHRhKSkgKyAKICBnZW9tX3Ntb290aChtZXRob2QgPSAiZ2xtIikgKwogIHRoZW1lX2J3KCkgKwogIHRoZW1lKGxlZ2VuZC5wb3NpdGlvbj0iYm90dG9tIikgKwogIHlsYWIoIkNpdHkgUmFuZG9tIEVmZmVjdCBJbnRlcmNlcHQiKSArIHhsYWIoIlJlZ2lvbmFsICgxMDBrbSkgU1NNIikgKyBsYWJzKHRpdGxlID0gImVCaXJkIikgKyBndWlkZXMoc2l6ZT1ndWlkZV9sZWdlbmQodGl0bGU9IlJlZ2lvbmFsICg1MGttKSBFbGV2YXRpb24gRGVsdGEiKSkKCmdnc2F2ZSgiY2l0eV9lZmZlY3RfcmljaG5lc3NfX291dHB1dF9fbWVybGluLmpwZyIpCmBgYAoKSG93IG11Y2ggdmFyaWF0aW9uIGhhdmUgd2UgZXhwbGFpbmVkPwotLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KYGBge3J9Cm1lcmxpbl9jaXR5X2RhdGFfbmFtZWQkcmVzaWR1YWxzIDwtIHJlc2lkKG1lcmxpbi5maXQpCmdncGxvdChtZXJsaW5fY2l0eV9kYXRhX25hbWVkLCBhZXMoeSA9IHJlc3BvbnNlLCB4ID0gcmVzaWR1YWxzKSkgKyAKICBnZW9tX3BvaW50KCkgKwogIGdlb21fc21vb3RoKG1ldGhvZCA9ICJsbSIsIHNlID0gRikgKwogIGxhYnModGl0bGUgPSAiTWVybGluIiwgc3VidGl0bGUgPSBwYXN0ZSgiQ29ycmVsYXRpb24iLCBjb3IobWVybGluX2NpdHlfZGF0YV9uYW1lZCRyZXNpZHVhbHMsIG1lcmxpbl9jaXR5X2RhdGFfbmFtZWQkcmVzcG9uc2UpKSkgKwogIHRoZW1lX2J3KCkKYGBgCgoKLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQpDaGVjayBBSUMKLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQoKYGBge3J9CkFJQygKICBnbG0oZGF0YSA9IG1lcmxpbl9jaXR5X2RhdGFfZml4ZWQsIGZvcm11bGEgPSByZXNwb25zZSB+IHJlZ2lvbl8xMDBrbV9zc20gKyByZWdpb25fNTBrbV9zc20gKyByZWdpb25fNTBrbV9lbGV2YXRpb25fZGVsdGEgKyBiaW9tZV9uYW1lICsgcG9wdWxhdGlvbl9ncm93dGgpLAogIGdsbShkYXRhID0gbWVybGluX2NpdHlfZGF0YV9maXhlZCwgZm9ybXVsYSA9IHJlc3BvbnNlIH4gcmVnaW9uXzEwMGttX3NzbSArIHJlZ2lvbl81MGttX2VsZXZhdGlvbl9kZWx0YSkKKQpgYGAKCmBgYHtyfQpBSUMoCiAgZ2xtKGRhdGEgPSBiaXJkbGlmZV9jaXR5X2RhdGFfZml4ZWQsIGZvcm11bGEgPSByZXNwb25zZSB+IHJlZ2lvbl8xMDBrbV9zc20gKyByZWdpb25fNTBrbV9zc20gKyByZWdpb25fNTBrbV9lbGV2YXRpb25fZGVsdGEgKyBiaW9tZV9uYW1lICsgcG9wdWxhdGlvbl9ncm93dGgpLAogIGdsbShkYXRhID0gYmlyZGxpZmVfY2l0eV9kYXRhX2ZpeGVkLCBmb3JtdWxhID0gcmVzcG9uc2UgfiByZWdpb25fNTBrbV9zc20pCikKYGBgCgoKCi0tLS0tLS0tLS0tLS0tLS0tLS0tLQpMREcKLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQpgYGB7cn0KYmlyZGxpZmVfZGF0YV93aXRoX2xhdCA9IGZldGNoX2NpdHlfZGF0YV9mb3IoJ2JpcmRsaWZlJywgaW5jbHVkZV9jaXR5X25hbWUgPSBUKQoKbW9kZWwyIDwtIGdsbShmb3JtdWxhID0gcmVzcG9uc2UgfiBJKGNpdHlfY2VudHJlX2xhdGl0dWRlXjIpLCBkYXRhID0gYmlyZGxpZmVfZGF0YV93aXRoX2xhdCkKZGF0MiA8LSBwcmVkaWN0KG1vZGVsMiwgc2UuZml0PVQpCnN1bW1hcnkobW9kZWwyKQoKd2l0aChzdW1tYXJ5KG1vZGVsMiksIDEgLSBkZXZpYW5jZS9udWxsLmRldmlhbmNlKQoKb3V0c2lkZV9zZSA8LSBiaXJkbGlmZV9kYXRhX3dpdGhfbGF0W2JpcmRsaWZlX2RhdGFfd2l0aF9sYXQkcmVzcG9uc2UgPCBkYXQyJGZpdCAtIDE1KiBkYXQyJHNlLmZpdCB8IGJpcmRsaWZlX2RhdGFfd2l0aF9sYXQkcmVzcG9uc2UgPiBkYXQyJGZpdCArIDE1ICogZGF0MiRzZS5maXQsXQoKZ2dwbG90KGJpcmRsaWZlX2RhdGFfd2l0aF9sYXQsIGFlcyh4ID0gY2l0eV9jZW50cmVfbGF0aXR1ZGUsIHkgPSByZXNwb25zZSkpICsgCiAgZ2VvbV9wb2ludChzaXplPTEpICsgCiAgZ2VvbV9zbW9vdGgobWV0aG9kID0gImdsbSIsIGZvcm11bGEgPSB5IH4gSSh4XjIpKSArCiAgZ2VvbV90ZXh0KGFlcyhsYWJlbCA9IG5hbWUpLCBkYXRhID0gb3V0c2lkZV9zZSwgc2l6ZSA9IDMsIHBvc2l0aW9uID0gImRvZGdlIiwgdmp1c3QgPSAiaW53YXJkIiwgaGp1c3QgPSAiaW53YXJkIiwgY29sb3IgPSAicmVkIiwgYW5nbGU9LTE1KSArCiAgZ2VvbV9wb2ludChkYXRhID0gb3V0c2lkZV9zZSwgY29sb3IgPSAicmVkIikgKwogIHRoZW1lX2J3KCkgKwogIHlsYWIoIkNpdHkgUmFuZG9tIEVmZmVjdCBJbnRlcmNlcHQiKSArIHhsYWIoIkxhdGl0dWRlIChvZiBjaXR5IGNlbnRyZSkiKSArIGxhYnModGl0bGUgPSAiQmlyZGxpZmUiKQoKZ2dzYXZlKCdjaXR5X2VmZmVjdF9yaWNobmVzc19fb3V0cHV0X19sZGdfYmlyZGxpZmUuanBnJykKYGBgCgoKYGBge3J9Cm1lcmxpbl9kYXRhX3dpdGhfbGF0ID0gZmV0Y2hfY2l0eV9kYXRhX2ZvcignbWVybGluJywgaW5jbHVkZV9jaXR5X25hbWUgPSBUKQoKbW9kZWwyIDwtIGdsbShmb3JtdWxhID0gcmVzcG9uc2UgfiBJKGNpdHlfY2VudHJlX2xhdGl0dWRlXjIpLCBkYXRhID0gbWVybGluX2RhdGFfd2l0aF9sYXQpCmRhdDIgPC0gcHJlZGljdChtb2RlbDIsIHNlLmZpdD1UKQpzdW1tYXJ5KG1vZGVsMikKCndpdGgoc3VtbWFyeShtb2RlbDIpLCAxIC0gZGV2aWFuY2UvbnVsbC5kZXZpYW5jZSkKCm91dHNpZGVfc2UgPC0gbWVybGluX2RhdGFfd2l0aF9sYXRbbWVybGluX2RhdGFfd2l0aF9sYXQkcmVzcG9uc2UgPCBkYXQyJGZpdCAtIDE1KiBkYXQyJHNlLmZpdCB8IG1lcmxpbl9kYXRhX3dpdGhfbGF0JHJlc3BvbnNlID4gZGF0MiRmaXQgKyAxNSAqIGRhdDIkc2UuZml0LF0KCmdncGxvdChtZXJsaW5fZGF0YV93aXRoX2xhdCwgYWVzKHggPSBjaXR5X2NlbnRyZV9sYXRpdHVkZSwgeSA9IHJlc3BvbnNlKSkgKyAKICBnZW9tX3BvaW50KHNpemU9MSkgKyAKICBnZW9tX3Ntb290aChtZXRob2QgPSAiZ2xtIiwgZm9ybXVsYSA9IHkgfiBJKHheMikpICsKICBnZW9tX3RleHQoYWVzKGxhYmVsID0gbmFtZSksIGRhdGEgPSBvdXRzaWRlX3NlLCBzaXplID0gMywgcG9zaXRpb24gPSAiZG9kZ2UiLCB2anVzdCA9ICJpbndhcmQiLCBoanVzdCA9ICJpbndhcmQiLCBjb2xvciA9ICJyZWQiLCBhbmdsZT0tMTUpICsKICBnZW9tX3BvaW50KGRhdGEgPSBvdXRzaWRlX3NlLCBjb2xvciA9ICJyZWQiKSArCiAgdGhlbWVfYncoKSArCiAgeWxhYigiQ2l0eSBSYW5kb20gRWZmZWN0IEludGVyY2VwdCIpICsgeGxhYigiTGF0aXR1ZGUgKG9mIGNpdHkgY2VudHJlKSIpICsgbGFicyh0aXRsZSA9ICJlQmlyZCIpCgpnZ3NhdmUoJ2NpdHlfZWZmZWN0X3JpY2huZXNzX19vdXRwdXRfX2xkZ19tZXJsaW4uanBnJykKYGBg